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ANALYSIS OF OPEN-LOOP SYSTEMS 
WITH PULSE-AMPLITUDE MODULATION 


Ya. Z. Tsypkin 


Systems with Class I pulse-amplitude modulation are examined. It is 
shown that in studying transient and steady-state processes in such systems 
it is possible to use the well-known results of the investigations of conven- 
tional pulse systems with Class II pulse-amplitude modulation obtained on 
the basis of the discrete Laplace transform 


An open-loop system with Class I pulse-amplitude modulation (PAM) 
can be represented in the form of a series-connected switch and a continu- 


ous circuit (Fig. 1) 


Fg, X@|  K@ AE 
FF (q,2) K*(qe—2) | 2° 2) 


Figure 1 


The key is closed periodically, at a given time, in such a manner that 
a series of pulses act on the continuous circuit with their peaks changing in 


‘accordance with the shape of the input signal during the interval of pulse 
duration (Fig. 2). 


a) 


Figure 2 


In investigating and calculating similar systems in the works of Far- 
manfarma [1] and Tou [2], special transforms were introduced (P—trans- 
form, T—transform)! which are essentially Laplace transforms of a mod- 
ulated series of pulses (Fig. 2b). 

In the present article it is shown that the study of systems with class I 
PAM or, as they are otherwise designated, systems with a switch, can be 
reduced to the study of conventional open-loop pulse systems — systems 
with Class II PAM, the theory of which has been amply developed [3, 4]. 
The mathematical method for investigating an open-loop pulse system — 
the discrete Laplace transform — makes it possible to carry out this in- 
vestigation by the most simple means. 

Let us designate the operating period of the key by T, and the duration 
of the closing of the key by yT. The output value of the system will be ex- 
pressed as 


Z(p)=K,(p)X(p), (1) 
where X(p) — is the expression for the modulated series of pulses, and 


K,(p) — is the transfer function of a continuous circuit. In relative values 
this equation is written in the form 


Z (q) = K (q) X (q). (2) 


Here, q = pT and 


ac 2 soi 
WE if K,(+). (3) 


Let us find the expression X(q). By definition 


2 (4) 


where 


<n+1 (5) 
and t = t/T. 
Equation (4), after substituting x(t) from (5), can be rewritten as 


co n+y 


X= f et faz (6) 


n=O n 
or, after substituting the variable t = n+ ¢ and changing the order of the 
summation and integration 


Xe liue “fatele “de, (7) 
0 


‘See Jury's book, Chapter 9 (GiaeS) rs 


But f[n + q = f[n, €] is a matrix function, and 


°° 


ve flac] —F* (ge) 
n=0 (8) 


determines the discrete transformation of this matrix function [4]. Thus, 
the onion for a modulated series of pulses may be represented in the 
form 


(9) 


Consequently, from Equation (2) and considering Equation (9), we ob- 
tain the equation for the investigated system (1) in terms of the general 
Laplace transform. 
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Z(q)=K(g) { F*(qele “ds. 
(7) (a) | (q,e)e ais 


The use of this equation requires operations on transforms which are 
functions of q and e4. In many cases this presents a complicated problem 
and special tables have to be prepared [1, 2]. 

In order to avoid these difficulties and to show that these systems can 
be examined on the basis of the theory of conventional pulse systems [4], 
let us proceed from a general Laplace transform to a discrete one. Ref- 
erence [4] shows that 


3 es) = ee \ a (q+2njm)é , 
Z (7,2) =D {Z(q)}= a Z(q+2njm) e a1) 


r=—o 


Using the discrete Laplace transform in Eq. (10), we obtain 


i nN 
Z*(q,2)=D{K(q) f Fee” ai i 
o (12) 


By changing the order of integration and using a D — transform, we 
get 


1 


aot yemex| “RT (an) DAKY(q)elLeyds 
Z* (4,8) J (0) D{K ge“ } a 


21t can be shown that this expression determines the P and T transforms 
(used in References [1 and 2]) using the discrete Laplace transiorm, 


Bn oe ie Eqs 
But in accordance with the theorem of multiplication by e an 


:. K (g,2 —+) we 
1) =| ee’ .K(q,l+e—4) leave (14) 


K* (g, 2) =D{K (a)}- (15) 


Consequently, the final equation of an open system with Class I 
pulse-amplitude modulation expressed in terms of a discrete Laplace 
transform is written 


if 
2* (g,2) = | Kx(q,e —») F* (q,4) dd 
0 


(16) 
or, in more detail, taking into account expression (14), 
6 5 ct 
Z* (g,e)= [ K* (ge—2) F*¥(g,Naddte ’ [ K*(g,1 +e—d F*(q,N ad 
0 = 
(17) 
at e<y, 
; 
and Z* (g,2)= | K*(g,¢—2) F* (9,0) dd 
! q ) F* (g, 4) (18) 


atl2eovy. 

Thus, knowing F(q) and K(q), from formulas of the type of Equation 
(8) or (11), we may obtain F* (q, €) and K* (q, €). Substituting these terms 
into Equations (17) and (18), we obtain an expression for the output value 
in terms of a discrete Laplace transform. 

Let us examine some special cases of the equations obtained for a 
Class I PAM system. 

1. For sufficiently small values of y (Fig. 3a) when the integrand 
can be considered constant, Equation (18) yields 


z* (9, e) = 7 K* (9,.£) F* (9,0). (19) 


Equation (19) is the equation for a simple open-loop pulse system. 
This result should have been expected since for sufficiently small values 
of y the shape of the pulses can be disregarded, and they may be con- 
sidered instantaneous with a sufficient degree of accuracy. 

2. At y = 1 (Figure 3b) 


*See Reference (4), p 184, where the theorem of multiplication by ea 
is given. 


1 
Z* (q, £) =) K* (g,¢ —) F* (g,4) dd. (20) 


xX at low 
values of C 


Lon a} 


Figure 3 


By using in this equation the inverse D - transformation [4], 


1 
—qe 


=] 
D'{2*(g,2)}= f Z*@e)e “de =Z() 
0 (21) 
we obtain an equation for a continuous system 
Z(g) =K (g) F (9). (22) 


It can therefore be concluded that for a system with Class I PAM the 
limiting cases will be an impulse system and a continuous system. 

The equation for an open-loop system with Class I PAM can be inter- 
preted in the following manner. The integrand 


Zo(qe —») = Kolqe — *) X* 9) (23) 
represents an equation for an open-loop, conventional pulse system where 
the pulse element reacts to the discrete values of the input value at the in- 


stant €=n+A [4]. 
Thus, in accordance with Equation (16) 


1 
Z*(@e)= f 2 (2— 4» os 


Inverting Eq. (24)*, we obtain 


z[n,<¢]= f Z)[n,e —A] a4, 


where 


4The parameter -integration theorem was used here. See Reto late. 230: 
Theorem 10. 


= wt 2 [me —2] when « >A 
Sie za zo[n—1, 1 +e —2] whene <h fe: 
and 
re 
29 [a3] = =f | Zi (9, 2) e" aq. (27) 


c—jr 
The method for calculating Equation (27) is covered in detail in refer- 
ence [4]. 
In more detail, taking Equation (26) into account, Eq. (25) can be re- 
written as two equations: 


Z[a—1,1+¢e—A]di (28) 


AbpCucer yin 


ae 
i 


z[n,2]= | 2)[n,2—] a), (29) 
f 
Be Deby 


Thus, in order to determine the equations and processes of an open 
system with Class I PAM, it is necessary to integrate Equation (27) from 
0 to y,, for conventional open-loop systems after substitution of « - A for 
Astor.0. 

This deviation makes it possible to use the well-known results ob- 
tained for open-loop pulse systems for the analysis of open-loop systems 
with Class I PAM. 

As an example, let us examine a system with Class I PAM whose 
transfer function for a continuous circuit is equal to 


k 


HONS lone as 


(30) 


Let us assume that an exponential effect is applied to the input of the 
system 


rc] 


fi(ij=e™. (31) 


The expression for this effect in terms of a general Laplace transform 
is equal to 


PA Ne rere (32) 


Proceeding to the dimensionless quantity q = pT, we obtain: 


K(q) = K( 5) =. tae (33) 


iP 1p qt 
; 1 1 
F(q)=—F,(4)= 34 
(q) Tr ea.) Ace (34) 
were a =a, 7; 6=8,T. 
By these expressions we obtain: 
1 ke ak: 
Ke 48) == { | == 
7°) gas) aa ete (35) 
eis —/D) [ l | = Ss a fe 
be) (Gn ete beth (36) 


Accordingly, the equation for a conventional pulse system (23) can now 
be rewritten: 


e27 — ae (a—3) A» 


R= —— 
(fe) (e7 —e F) e (37) 
Reverting to the transforms from the inversion, we obtain 


k —a (n+l) —B (n+1) —1e (a—B)A 
—h = ———, (e —e e e 

Zo[n,€ ] a ( ) (38) 

Substituting this value into Equations (28) and (29), after integration 


and elementary transformation we finally obtain 


— 32 


z[n,e] =k sacra elle [(es a) ae 
(e* — ec?) (2—8) 
or [(e* 2 es) ) aC, (eo Ms ate (39) 
(€< 7) 
and 
(a—P)t _y —a(n+1) —A(n+1)) ,— ae 
z[n,e] =k : e —e Merny 
(et — e~* (a —8) (20) 
(e > 7) 
When y is small, from Equation (40) we obtain 
Cin ce —a (n+1) —B (n+1) 
Z|n, ej = SME ace, ) 
[n, ©] = ky aaapee ( 8) 
which corresponds to the process in an impulse system. 
When y = 1, from Equation (39) we obtain 
1 —a(n+e — B(m+e) 9 
z[n,e]=k =| (e Catia ee eee nan): (42) 


which, atn+ € = ‘t, corresponds to the process in a continuous system. 
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NOISE PROOFNESS IN INTEGRAL RECEPTION 
OF SIGNALS UNDER THE INFLUENCE 


OF FLUCTUATING AND SINUSOIDAL NOISE 


N.L. Teplov 


A comparison is carried out of noise proofness of integral and con- 
ventional narrow-band reception of pulse signals under the influence of 
fluctuating and sinusoidal noise+. 


INTRODUCTION 


Asa rule, methods of integral reception are compared with methods of 
conventional narrow-band reception. Such a comparison is a natural one 
because integration (averaging), as a physical process, is equivalent to 
filtration (smoothing out). 

Narrow-band reception is usually understood to mean the reception of 
signals with the most adequate passbands which correspond to maximum 


‘The evaluation of noise proofness in the indicated methods of reception for 
other forms of noise is given in a separate article. 


noise proofness of reception. The values of the indicated bands are not the 
same for different forms of noise. For example, in the case of fluctuating 
noise, the most suitable bands for the reception of single pulses are the 
so-called optimum frequency bands [1]. 

A characteristic feature of conventional narrow-band reception (re- 
ferred to hereafter simple as narrow-band reception) is that the filtering 
of the signal and noise is accomplished by means of filters with fixed con- 
nections. Accordingly, the selectivity of the latter is determined by their 
static resonant characteristics and passbands. 

The most suitable values for the filter bands, as indicated above, cor- 
respond to the conditions (and consequently, the most suitable values for 
the filter passbands) for identical forms of noise are found to be different 
during the reception of single pulses (following one after another with 
large gaps in between, for example, as in radar) as compared with the re- 
ception of a continuous series of pulses (for example, in radio telegraphy). 
In the first case, mutual overlap of the pulses due to transients in the fil- 
ters practically never takes place. Whereas, in the reception of a contin- 
uous series of pulses, the appearance of each pulse takes place in the pres- 
ence of residual oscillations in the filters fr»m preceding pulses. As a 
rule, the most suitable frequency bands in the latter case are found to be 
wider ones: they must, as a compromise, satisfy two contradictory re- 
quirements — the assurance of the highest possible increase in signal over 
noise in the filter, and minimum interaction between adjacent pulses. 

Thus, for conventional narrow-band reception the conditions for re- 
ception with minimum probability of error are found to be different (do not 
coincide) for different forms of noise, and depend also on the pulse-repiti- 
tion period in the total sequence. 

Let us now examine the features of integral signal reception. 

It is understood that by integral reception we mean reception with ac- 
cumulation of the signal and noise only during the interval of signal dura- 
tion. Integration in this sense can, in practice, be realized in filters (cir- 
cuits) with very narrow bands, 2 switched on and off in phase with the rise 
and fall of each signal pulse. 3 Switching filters (circuits) with very- 
narrow frequency bands shall hereafter be referred to as integrating filters. 

We may formulate the following distinctive features of the integral re- 
ception method: 

: 1. The effect of noise on the integrating circuit takes place only in the 
- interval of signal duration (in contrast to reception on narrow-band filters 
with fixed connections in which the voltages and currents are generated 
continuously by the noise during its full operating period). 


2 By very narrow-band filters we mean those in which the time for the 
establishment of steady-state conditions is known to be greater than 
the duration of a single signal pulse. 


3% radiotelegraphic apparatus known at this time, the switching on and 
off of the integrating filters in phase with the limits of the unit pulses: 
is achieved by means of forced quenching of the residual oscillations in 
the filters at the end of each pulse; in this way, null conditions are 
established in the integrating circuits at the beginning of each fol- 


lowing pulse. 


2. During integral reception, the mutual superimposition of adjacent 
pulses, caused by the long duration of transients in the narrow-band fil- 
ters, is eliminated; this is achieved through artificial (forced) quenching 
of the residual oscillations in the integrating circuits at the end of each 
pulse separately. 

The above-indicated features of integral reception lead to a new quan- 
titative correlation of signal and noise which is different from the analagous 
correlation for conventional narrow-band reception. The qualitative dif- 
ference of integral signal reception from conventional narrow-band recep- 
tion stems from the fact that its frequency selectivity is determined not by 
the static characteristics as in narrow-band reception, but by the dynamic 
resonant characteristics and bandwidth of the integrating circuits. 

The overall advantage of integral reception methods is their greater 
noise proofness. At the same time, the conditions for highest noise proof- 
ness for various types of noise are satisfied simultaneously (mutually 
compatible); whereas with conventional narrow-band reception, the maxi- 
mum attenuation of different types of noise is attained without identical 
conditions for the various noise signals (especially with different values 
for the filter bandwidths. ) ; 

We shall also note that noise proofness of integral signal reception 
does not depend on the pulse repetition rate of the sequency (since with in- 
tegral reception the interaction between pulses is eliminated). 

The noise proofness of integral reception of pulse signals is evaluated 
below, with integration of signal and noise in the radio frequency channel 
of the receiver (to the detector stage). 

Figure 1 shows a block diagram of the receiver with integration to the 
detector; the integrating circuit is designated as the integrator for sim- 
plicity. 

A corresponding block diagram for a narrow-band receiver is given 
in Figure 2. 


Common 
receiver 
channel 


Input Detector 


Figure 1 


Common 
receiver 
channel 


(Go) 


Narrow- 
band 
filter 
(Wy) 


Input Detector 


Figure 2 


The whole radio frequency channel (up to the detector), which also in- 
cludes the integrator or the narrow-band filter, is tuned to the signal fre- 
quency Ww, = Wo. 


0 
In the general case, we shall represent the signal in the form of a 
series of radio-frequency pulses with rectangular wave-shape 


A(t) =asinat when0<t<r (1) 
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where a and T are the amplitude and duration of the unit pulses. 

It is further assumed that the signal is not noticeably distorted in the 
common (relatively wide bandwidth) channel of the receiver. 

The pulse repetition rate of the series may be varied. 


NOISE PROOFNESS OF INTEGRAL AND NARROW-BAND 
RECEPTION UNDER THE INFLUENCE OF FLUCTUATING NOISE 


Integral Reception. Fluctuating noise at the integrator input is repre- 
sented by the oscillation 


B(t)= 6 (é)sin [oy 4+ 6 (], (2) 


where b(t) and @(t) are the random amplitude and phase. 
The quadratic for the effective value (intensity) of the fluctuating 
noise (2) is equal to 


[BOP=o? =A fore, (3) 


where Afp¢r is the effective frequency band of the receiver channel at the 
integrator input; v9 — is the noise intensity in a 1-cps band (specific noise 
intensity). 

The oscillation (2) may be divided into the components 


B(t) =x () sino, t+ y (0) cos af, (4) 
where 
x (t) =6(t) cos 6 (A) 


yY=6 /sin6 (4) (5) 


The probabilities of the amplitudes x(t) and y(t) are normally distribu- 
ted with an average value equal to zero, 


x=yO=0 iS) 
and dispersion equal to 
[OP =) OF =? = Sete (7) 
The oscillation (4) is also rep- y 
resented by the vector diagram in 
Figure 3. 
Let us determine the intensity % 


of the fluctuating noise at the inte- 
grator output. 

In accordance with the expres- 
sion (4), the fluctuating noise at the By 5 
integrator output is a sum obtained 
by integrating the oscillations of 
amplitudes x(t) and y(t). The am- Lary ie ' 
plitudes of these oscillations after their integration in the interval of signal 


duration T is equal to: 


Figure 3 


Lal! 


x= [ x(pat 
0 


y= fy ({) dt 


0 


eS 


(8) 


The resulting amplitude of the fluctuating noise at the integrator out- 
put is equal to the geometric sum of the amplitudes (8) [See formula (4) 
and Figure 3]. 


bhe= p JO aes, (9) 


The amplitudes of the oscillations x, and yr are random values with 
normal distribution (as are the instantaneous values of the amplitudes 
x and y). The intensity (dispersion) of the fluctuations of the amplitudes 
X7 and y; are shown to be equal to [3] 


= = 9 


2 yp an x(i) dt)’ = fi y (Nai)! =o =a ans 


The mean quadratic value of the amplitude of the fluctuating noise at 
the integrator output is determined in accordance with Formulas (9) and 
(10) 


bX Te oy (11) 


The dispersion of the fluctuating noise (the quadratic of its effective 
value) at the integrator output is equal to 


o 
=e = ve, 
[BOk= 5 0 (12) 


The amplitude of the signal at the integrator output is equal to 


a, =| adit=az. 
J (13) 


Let us take the power ratio between the signal and noise at the inte- 
grator output [which, quantitatively, is equal to the ratio of the quadratics 
of the signal (13) and noise (11) amplitudes] 


The power ratio between signal and noise is usually called the signal- 
to-noise ratio. 


Accordingly, the signal-to-noise ratio at the integrator output is equal 


to the relationship between the signal energy ( “— ) and the specific 
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noise intensity (v ah and does not depend on the frequency band of the re- 
ceiver channel up to the integrator. 

The signal-to-noise (14) corresponds to the maximum noise proofness 
of signal pulse reception. On the strength of this ratio, (14) may be con- 
sidered the limiting value. 

Thus, integral reception assures maximum noise proofness in pulse 
signal reception under the influence of fluctuating noise. 

Narrow-Band Reception. The maximum signal-to-noise ratio for fluc- 
tuating noise in conventional narrow-band reception is realized with so- 
called optimal frequency bands of the filters. 

Depending on the shape of the resonance curve of the filter, the values 
for optimal frequency bands under the action of single pulses are found 
within the limits 


0.7—1.4 
Af eft opt i = , (15) 


where T is the pulse duration of a signal of rectangular waveform. 
The optimal frequency bands correspond to the following maximum 
signal-to-noise ratio 


2 = ¢ a2t 
Hage = 0.82 ats 


The relationship (16) holds true for any filter with an optimal fre- 
quency band. This position is illustrated below by examples corresponding 
to two extreme cases: 1) the case of a single oscillatory circuit, and 
2) the case of an ideal filter. ; 

The amplitude of a signal at the output ofa single resonant circuit (src) 
with radio pulses of rectangular waveform acting on its input is established 
according to the law 


@our= (1 — eet 4), 


where Afarr is the effective frequency band, related to the circuit attenu- 
ation factor by Afyrp = so the transmission factor of the circuit at res- 


onant frequency is taken conditionally as equal to unity. 
The maximum value of signal amplitude at the circuit output, at the 
instant t=7T, is 


—2Af t 
A utmax — 2 (dee ei iF 


The signal-to-noise ratio at the circuit output when t = T is equal to 


Milena ett SP OM ee cs a 


Qe M eff 2 Meffi* (17) 


hi = 


stc 


The maximum value of the ratio (17) corresponds to the optimal fre- 
quency band for the circuit which is determined from the conditions: 
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; ae a | 0 
d (AS eff) Meff * “Veff ~ “Seffopt 
By taking the derivative, we obtain the equation 
—24f eff opt * A Wea) 
ie be S eff opt *) , 


the solution of which gives 


0.65 
Af. ot 
eff opt (18) 
Substituting the value for the optimal frequency band in the expres- 
sion (17), we have 
9 FS Xe 
hopesre = 0.815 ——. (19) 


= 


The maximum value of signal amplitude at the output of an ideal filter 
is equal to [1] 


») 
Goatees aie Si ia Af), 


where Af is the filter passband: Si is the symbol for the integral sine. 
The signal-to-noise ratio is equal to 


bes (Gaf 


“aif (20) 


The optimum frequency band for an ideal filter will have the value [1] 


1.37 
Afopt = aa 


The substitution of this value for the frequency band in the formula for 
the ratio (20) gives 


2 of nendes 
opt F == 0.825 >? 


2vp (21) 


h 


which practically coincides with the formula (19) for a single oscillatory 
CIV CUM, 

The relationship of signal-to-noise ratios for integral reception and 
narrow-band reception with optimal frequency bands is equal to 


h, \ SSATiC 
Geet. = 1.22. a 
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Thus, during the reception of 
single pulses (followed by long 
gaps), integral reception and con- 
ventional narrow-band reception 
with optimal frequency bands are 
practically equivalent. 

The advantages of integral 
reception show up more effec- 
tively for a continuous sequence 
of signal pulses (such as take 
place in radiotelegraphy). 

Conventional narrow-band reception with optimal filter frequency 
bands is characterized by the relatively high duration of transient phe- 
nomena. This leads to the superimposition of residual oscillations from 
the preceding pulse on the incoming pulse (Fig. 4). The indication of each 
incoming pulse now results not from the absolute value of its level but from 
the relative difference in level Aa. The optimal frequency band in the case 
under examination is determined by the conditions for the highest ratio be- 
tween the difference in level of adjacent pulses and the fluctuating noise 
intensity. 

For example, in the reception of frequency-modulation signals, Aa is 
the difference voltage in off-and-on keying of messages at the output of the 
frequency detector (discriminator). 

Let us determine the optimal frequency band, taking into consideration 
the superimposition of pulses, for a single oscillatory circuit 


—2/ off = —23foff = ae cae = ~ 


—a(l—e 


== (71 (5) Be peel eit “ye 


ws Ss 


The signal-to-noise ratio in the case under examination will be de- 
termined from the relationship 


Pee Se 2 ea, (Nees cone ©) 
22 AJ/eff 22 Afeff (23) 


The maximum signal-to-noise ratio (23) corresponds to the optimal 
frequency band 


Af ett opt Spee (24) 


Substituting (24) in Formula (23) we obtain 


Ae eae 
opt 22 (25) 


It can be proved that in receiving a series of pulses the optimal 
frequency bands will be approximately twice as wide as the optimal 
frequency bands for single pulses, and the maximum signal-to-noise ratio 


will be equal to 


15 


| CO) ; 
or 23 (26) 


Taking the relationship between (14) and (26), we obtain 
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a, llega (27) 


Thus, with integral reception of a continuous series of pulses, the 
gain in signal-to-noise ratio in comparison with conventional narrow-band 
reception is shown to be practically double. 

The verifications carried out above relate to the reception of a con- 
tinuous series of pulses of equal amplitude. The gain in noise stability for 
integral reception was shown to be particularly significant when receiving 
signal pulses with amplitude variations (in particular, signals subjected to 
fading). Under these conditions the optimal frequency bands will be even 
wider and signal-to-noise ratio, smaller. 


NOISE PROOFNESS OF INTEGRAL AND NARROW-BAND RECEPTION 
METHODS UNDER THE INFLUENCE OF SINUSOIDAL NOISE 


Integral Reception. We shall write sinusoidal noise in the form of an 
oscillation 


En =e, Sin, t= e, Sin (w, + Ao), (28) 
which is detuned relative to the signal frequency by Aw = uy, — Wo. 

The duration of the sinusoidal noise T,, can be practically any value, 
but it is necessary to differentiate between the cases: 


a) sinusoidal noise of long duration exceeding the signal duration, 


heen 
b) short-time sinusoidal noise, when 
i Pal 


We shall evaluate here the effects of sinusoidal noise of long duration 
(‘non-attenuated" sinusoidal noise). This permits the expression of the 
frequency-selectivity properties of the integrator in a universally- 
adopted sense. 

Let us rewrite the expression for sinusoidal noise (28) in the form 


Ey = en CoS Awlsin wt + en sinAwfcoso,t—o (Ao, f) sin w, ¢ + 


“bE b (Aw, t) COS Ot for O < t < “pn rl (29) 


The oscillation (29) corresponds to the vector diagram in Figure 5. 
The amplitude of the sinusoidal noise at the integrator output is 
equal to 
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ep (Ao). = ae t) dt)’ a [9 (ao, t) dt’ = 
0 0 


=é, V [feo Aot at)" oa f sin Aot dt if 2 


— ou / — [sin? Aw: + (1 —cosdo:)?] = 


2 4 eho) 
2=G = (i = eos Ate) = meiner 
a Vy Aw? ( ) en hee 2 fi 
Accordingly, 
Aw 
Sint sin ai © 
en (40), =e ain = 
w 
By, Fave: 


For Aw = 0 (at signal frequency) 
en (0), =+. (31) 


Thus, the transmission factor for an ideal integrator at resonant fre- 
quency is equal to 


h@=s 


The relationships (30) and (31) determine the equations for the res- 
onance characteristic of the integrator 


sin 22 
ep (Ae), ! 2 
OR Se emer ear (ue 
: (32) 
where 
en (YY), sin Aft 
Ae ee 
y( d); €,(0),. Ast (33) 
_&® 
where Af = a ae 


The resonance characteristic of the integrator, calculated by Formula 


(33), is shown in Figure 6. 
Let us determine the value for the effective frequency band of the in~ 


tegrator 
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Figure 5 Figure 6 


Accordingly, the effective frequency band of the integrator is equal to 
the inverse value of the integration time (equal to signal duration). 

Narrow-Band Reception of Signals. In evaluating the frequency se- 
lectivity of narrow-band reception, it is necessary, above all, to keep in 
mind the limitation of signal duration. Because of this, with a decrease 
in the filter passband, we shall have, together with the limitation of noise 
effects, a simultaneous decrease in signal amplitude at the output. A 
proper criterion for evaluating the actual selectivity of the frequency- 
selective channel (of the filter) in the case under examination can be the 
relationship, at its output, between the maximum amplitude of a signal 
of given duration and the amplitude of the sinusoidal noise in a steady-state 
condition (as in integral reception, we have examined here the effect of 
sinusoidal nvise of sufficiently large duration, greater than the time re- 
quired for steady-state conditions in the frequency-selective system). It 
is assumed that at each given moment only one sinusoidal noise will be 
effective at one of the frequencies (the case of separate sinusoidal noises). 

A complete evaluation of the effects of sinusoidal noise in each con- 
crete case can be carried out if we know all the statistical facts of the si- 
nusoidal noise (in particular, if we know the distribution of the sinusoidal 
noise by frequency and level). A detailed treatment of the problem of sta- 
tistical structure of sinusoidal noise (which can be quite different depend- 
ing on actual operating conditions of a given communications system) is 
beyond the limitations of the present article. Here, we shall have in mind 
sinusoidal noise with an equal-probability distribution in frequency and 
with practically any values of amplitude. In other words, noise at any fre- 
quency is assumed to be of equal probability and, on the average, none is 
predominant over the others in an energy relationship. 

Let us introduct the concept of signal-to-noise ratio at a given fre- 
quency 


where e! is the mean statistical value of the sinusoidal noise amplitude 
at the given frequency. 

In accordance with the above, the mean statistical value of the sinus- 
oidal noise amplitude is assumed to be the same at each frequency. 

Then, the average, with respect to the full range, power of the sinus- 
oidal noise acting on the receiver with an effective frequency band Af ff 
will be equal to (with a 1-ohm resistor) 


2 x9 
en 


en t 
"g a ee Afere = 5 Afett (35) 
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2 : 
Here & Nos is the quadratic for the effective value (for a suf- 


n 
2 
ficiently large period of time) of sinusoidal noise acting on the receiver at 

one of the frequencies. 

Let us calculate the signal-to-noise ratio separately for integral and 
narrow-band reception. 

For integral reception we obtain on the basis of Formulas (13), (35) 
and (30) 


(36) 


By comparing Formulas (35) and (3), it is easy to conclude that the 
determination of maximum excess of signal over the aggregate sinusoidal 
noise during narrow-band reception can be reduced to the problem, solved 
above, for sluctuating noise. 

Based on this fact, we may write: 


2 0.82 
Font = 0,82 par (37) 
2 ee 05 
Rog = 0.5 02’ (38) 
where oe and hg ot represent the signal-to-noise ratio for sinusoidal 


noise, which correspond to the optimal frequency bands for the reception 
of single pulses in a continuous series. 


Accordingly, 
h 2 
—*) =1,22 
( me (39) 
and 
he V2 2.0 
Gor get’ (40) 


We shall recall that the results obtained relate to the case where the 
pulses of the received signal have an equal amplitude. If the pulses in the 
received series have different amplitudes, then the noise proofness of 
conventional narrow-band reception, in comparison with the integral- 
reception method, is decreased to a greater degree. 
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CALCULATION OF NON-LINEAR DISTORTION 
IN WAVEGUIDE TRANSMISSION LINES 


Yu. A. Yerukhimovich 


An exact solution is given for the problem of any number of inhomo- 
geneities in a waveguide line using the method proposed in reference [1]. 
A comparison with the approximate solution given in [1] is made. Some 
results of the exact solution are presented and an example of practical 
calculation is given. 


In [1] as well as in previous works of the author, an approximate solu- 
tion was given for the problem of long feeders when the number of inhomo- 
geneities, M, is large. In the present article, an exact solution is given 
for any number M, within the limits of applicability of the method. 

We shall consider the results of reference [1] as well-known, and 
hereafter, for uniformity, we shall use the notation introduced therein. 

It should be noted that another method for the solution of this prob- 
lem was used in reference [2, 3]. However, the results obtained by this 
method cannot as yet be used in practice for sectionalized waveguides. In 
accordance with [1], we shall assume that the waveguide (feeder) has M 
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inhomogeneities distributed at equal distances from each other, and each 
one of which has the same reflection factor. 


CALCULATION OF HARMONIC NOISE IN WAVEGUIDE 
LINES FOR ANY NUMBER OF INHOMOGENEITIES 


The noise power at the qth harmonic at a frequency wa, relative to 
the fundamental, is given by the equation (3) in [1]. Under the conditions 
stated above, P, may be written as: 


q 
P,=2( caves bint g Oe et? >) ea 0) 
M—2 M—1 
(SMe, 
EGE 9 


(1) 


Here 6 = 2Aw = ; Wa is the audio angular frequency of the test- 


tone modulation, Aw — is the peak frequency deviation of the multi- 
channel signal; L = M/ — is the total length of the waveguide; 1 — is the 
length of one of its sections; Vp — is the group propagation velocity of the 
multi-channel signal in the waveguide; ro, rjq — are, respectively, the 
reflection factor moduli for the generator and for the load (antenna); r — 
is the reflection factor modulus for each irregularity; Iq(z) — is a Bessel 
function of the first kind and qth order. 

The summations are replaced in [1] by the corresponding integrals to 
which they convert with unlimited increase in M. It is assumed that for the 
final M, its value must be sufficiently large so that the error in the cal- 
culations will be small. However, in the work cited, it is not indicated 
which value of M would be considered acceptable. Neither does it take into 
account that from the summation in (1), it is possible to separate out 
terms proportional not only to M* or M, but also to lower degrees. A 
study of Formula (1) showed that the summations it contains may be cal- 


culated of any value of M. 
For intermediate calculations we shall use the following relationships 


[5] 
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the final term contains M or M2. 
In order to avoid performing duplicate calculations we shall consider 


attenuation in the feeders. In accordance with [1], the first summation in 
Equation (1), taking attenuation into account, may be represented as: 


M-1 


S,= > Pie Her 0). 


Ff fie 
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Here, @ — is the attenuation factor of a wave in a waveguide of unit 
length. S07" 

Let us expand the expression under the summation sign into a series 
and change the order of summation: 


co o M-—1 


(—42l)s (/_ 8 \V | pts 
on 2 ps Ds Aan s! or: ) sabal, 


=0) s=0) m=1 


and y(M — 1) = ——— =. 


6 
2M 
Calculation of the summations with respect to m gives 


Let us insert the symbols y = 


M-1 
pts nets Dp pert i pee 
Y : — ae 


i Basel 2 -. 12 


m=1 


Using the integral and differential forms of the function, we may 
write 


p+s 7 d 
Fi ! 


! pts 
= x 
- = 2 Wendx 


m=!) 


Let us substitute this relationship inS], and take advantage of the 
feasibility of revers ng the order of the integration and differentiation 
operations and the summation operations (in an absolutely-converging 
series). The summations appearing in the resulting expressions are easily 
calculated by the formulas derived above; as a result we find 


4.1 4al 
~ ——-* 


Si=— i etnlwod? (2x) dxictarse ohy He (2a)ych 
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By substituting x = 6 M + under the integral sign, defining 
Vig : : 
6 iA cai Tae 6 (Fig. 1) and using the average value theorem for inte- 


gration of the product of two functions (as indicated in [1]), we find 


1 1 
Al —4a 
Si=(M—1) fe dy (RO x)dx+ te" 20,) + 
0 0 
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Figure 1 


For comparison with the results in the work [1], we shall use the 
symbols employed there: 


1 
FLO=hO FO) = J he Ox) dx 


1 
F; (8) = J (1— x) (6x) dx, 
0 


-4aL 
Abo = -4aL = 1—e 
and also the quantities Ky e and Ky yaa The symbol 


Kim means that in place of the quantity L, in the expression corresponding 


to K; we must substitute the quantity L),. Finally, we shall obtain 


ss (M — 1) F, (0,,) Kya, + O50 F: (8, ) Aya + 


—4Aal 
eet) 
24M dx Ons 
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Now we shall find the value of the double summation in (1). 
After a simple transformation, S5 (taking into account the wave atten- 
uation) may be represented in the following form: 


Using the relationships deduced above, we find 
oo oo M—2 M—n-1 


Se eG a) 


6=0) 1 y=0' In=1 s=l 


Let us determine, first, the summation with respect to s 


M—na-1 
© pivtl pty (PEN Op actt 
Sa ee ane ee 
s=1 
hereu =M—n-—l1. 
Further, we shall find the values inside the summations. In calcu- 

lating the summations with respect to n, we shall use the following auxil- 
iary relationship: 
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Let us designate z= 7(M — 24 ae at 


see Figure l). Matas 2 

Considering values not greater than the second order of smallness 
relative to 1/M, we shall have, after determining the sums according to 
n, 


(For the value @ M : 
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In order to simplify this expression, we shall use the following rela- 
tionship: 
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After transforming the integrals to fixed limits and using the theorem 
on the average integral for separating out the dependence of So on the wave 
attenuation, we have 


Sy (M— 2° Fy (0m) Kam + (M—2)F, (On) Kam + Fy On) Kina + 


je eS Pe: 
wae (22) ae 
z <n OM (3) 
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Here, Ly = re a L (Fig. 1); the symbol «’» for | ieee means that in 


ihe place of L in the explicit form of this function we should substitute 
M: 
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Now we shall write the final expression for Pg: 
@q \2 
Po=2(q 22) {BA Fi (8) Ky + (M—1) 2 (B+) Fe Oy) Koy + 
( A 5 ' ' 
+ (M — 2) F,(6,,) Koy + ToL! (On eGrse 


1 ' ' 
+> OG +E On) Kin + (M— 2)? 4 Fs (8y) Kom |. Pi 


The functions Ky and K, are shown in Figure 2, in accordance with 
those given in [1]. 

The argument of the function F,,(@) in reality is found to be dependent 
on the number of sections in the waveguide. The same applies also to the 
effective length of the whole waveguide. 

An evaluation of the accuracy of the calculation performed by means 
of (4) may be obtained from the formula 
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Figure 2 


auiae quantity 4 Ba designates the expression within brackets in formula 


4). 
y For M—> ~, we have Ly L;: L yr: 6M—é@ and this formula is trans- 


formed into (4) from [1] with the following two comments: 

1) From Formula (8) in [1], it follows that in order to calculate the at- 
tenuation F3(8) in (4), it is necessary to multiply by Kg, which differs from 
Ky. However, this conclusion is in error. Actually, it follows from [1] 
Appendix 1) that the function F3(8), with attenuation taken into account, is 


determined by the expression 
—4dal (y—*x) 


11 
F,()=J fe P[8(y — x)] dx dy. 


By a formula derived in this same appendix, it is identically equal to 
: 4aLl 
fe 0/5 (6x) dx. 
0 


By applying to this integral the theorem on the mean integral from the 
product of two functions, we find that in calculating the attenuation F3(6), 
it is necessary to multiply by 


1 
(ede ne 
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2) From the results of a precise summation of S, and Sp, itis clear 
that in [1], during the substitution of integrals for the summations, a 
number of terms were not taken into account. This circumstance is the 
principal shortcoming of the method used in [1] for calculating the sums. 


NOISE IN NONLINEAR TRANSITIONS 


Of basic interest are the formulas which permit the computation of 
the relative or absolute value of noise in the nonlinear transitions during 
the transmission of multichannel messages through the transmission line. 

By using the results and data obtained in [1] and in previous works, it 
is easy to write the general formula analogous to (7) in [1]. 

For this purpose it is sufficient to note that @ m not equal to 6, the 


quantity 6 should be replaced by the value B yy = (— B and @ M by 


' Q' 
Bip = (—g—)?B, ive. 


Pu=(—-) 8 (6) 
n= (SY 8. (7) 


Now the general expression for the relative noise-power level in non- 


linear transitions, analagous to (7) in [1], for the highest frequency of the 
multichannel signal, may be written as: 
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where x — is the utilization factor of the total number of channels (x ~ 1/4). 
In accordance with [1] 
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B =472 x2, T= 


x — is the effective value of angular deviation of the multichannel signal, 


is the signal transmission time in the waveguide; 


which is equivalent to == for the test signal. 
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Although now the expression for G? has become more complex because 
of the addition of new terms, all of them are determined by the same func- 
tions F,,(8) as in work [1], but with new values of the argument dependent 
on M. The functions F,,(8) have been reproduced here in Fig. 3 from data 
in [1], and may be used directly for calculations by Formula (8). 


Figure 3 


The accuracy of the calculations carried out by means of Formula (8) 
can be determined from the expression 
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Here { G?} denotes the quantity within brackets of formula (8). Dif- 
ferentiating with respect to the parameter, we find 
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This function is shown in Figure 4. 
The basic result of the precise solution is worth noting. The dependence 


of the arguments of the functions F,(8) and Ky(@L) on the value of M leads 


27 


to an exaggerated evaluation for the degree of influence exerted by these 
functions on the overall noise level. A particularly significant effect of 
this factor will be observed in the region of rapid changes in F,,(8) of the 
function 8, i.e. when Tx~ 1. At the same time, as actual calculations 
will show, the value for noise determined from data in [1] may appear to 
be of an order higher than that obtained from precise calculations. For 
large values of Tx and M, the difference between formulas (8) and (7) in 
{1] should level out. 

From Figures land 2, it follows that for slightly larger values of M, the 
wave attenuation in the waveguide has a greater effect on the value for 
relative noise level than is indicated in [1]. ! 


Ox om 
In the case when M is so large that and ———— approach 1, 


and B does not lie in the region of rapid changes in F (8), the formula for 


noise of nonlinear transitions in the approximation which is used in the 
work [1] (M > 1) acquires the form 
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As an example of the application of the obtained results to practical 
practical purposes, we shall present the results of calculations of noise 
level at nonlinear transitions of a waveguide transmission line for a 
typical radio relay line. The magnitude of the noise is taken relative to a 
power of 1 mw at the mean audio frequency of the highest channel of the 
multichannel message, and is expressed in picowatts. The scale of the 
curves takes into account the psophometric (circuit noise) factor, for 
which the results of calculations by formula (83) [sic] are multiplied by 
the coefficient k® = 0.56 [2]. 

The waveguide consists of M sections each of which has a length 

= 4.5 m. The noise relationships were computed in the function M for 
various values of the parameters rg and ry, and for simplification, it 
was considered that the reflections from the waveguide terminations were 
equal to each other, i.e. rq = ryg. The magnitude of the reflection from 
an individual junction was considered to be an independent parameter 
Both rp and r are expressed in percentages. Figure 5 shows the results of 
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calculations for the following data: number of channels N = 240; attenuation 


aa db 
O= 2.5 joo mo X= 550 ke and the angular frequency of the highest chan- 
nel fy = 1052 ke. 
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Figure 5 


The numerical values for the parameters in Figure 6 are: N = 600; 
Sa Sue 

OO Maa? 
curves, the following rules may be established. With small reflection 
factors for the unctions, the curves have clearly indicated noise peaks 
located near the noise peak for a smooth (continuous) waveguide (r = 0). 
With an increase in r, the positions of the peaks are displaced in the direc- 
tion of longer waveguide lengths (larger M). 

For small values of reflection factor for the waveguide terminations, 
the principal role will then be played by the values of reflection factor for 
the junctions. For small r, the noise level for this calculation will be 
linear with increase in waveguide length. The greater the value for N, with 
a constant value for ro = r,,, the shaper will be the noise peaks. 

On the basis of fhe computed curves, it is possible to determine 
the conditions under which the noise level in a system will not exceed 
acceptable limits. 

1) The results attained have expanded the possibilities of using 
the method written up in [1] and previous works, and makes it pos- 
sible within its framework to analyze sectionalized transmission lines 
with any number of links. The errors in [1] have been eliminated. 

2) The deduced data makes it possible to make an evaluation of the 
accuracy of the calculations performed. 

3) Existing recommendations concerning the computation of noise re- 
sulting from nonlinear transitions (see, for example, [4]), require cor- 
rection. 


x = 880 ke and fy = 2540 kc. By analyzing the resulting 
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A SIMPLE ELECTRONIC INSTRUMENT 
FOR STATISTICAL RESEARCH IN THE FIELD 
OF RADIO WAVE PROPAGATION 


N.S. Stanulov 


A description is given of an electronic instrument for the direct de- 
termination of the average value of field intensity. An ordinary de electric 
meter was used as the indicator of average values. With known distribution 
of the instantaneous values of the electromagnetic field intensity, the in- 
strument may also be used for determining the corresponding median val- 
ues. 


In the theory of radio wave propagation, the study of ionospheric 
propagation occupies an important position. Under the conditions of this 
type of propagation the amplitude of the resulting field at the receiving 
point changes in a random manner, i.e. fading occurs. The reasons for 
these oscillations in field intensity are the continuous variations in the 
electron density of the ionospheric layers, oscillations in their altitude, 
and interaction between waves which are propagated along many different 
paths [1, 2]. Quantitatively, the fluctuations in field intensity may be de- 
fined by some random function of time E(t). Figure 1 shows a possible form 
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of this function (recorded in the course of a certain interval of time). 
Characteristic for this type of recording is the fact that each time the ceo 
cording apparatus is switched on, the function has a different form. This 
shows that the oscillations in field intensity represent one of several pos- 
sibilities E,(t) (kK=1, 2, 3, .. -) of the random function E(t). 

For a qualitative and quantitative evaluation of fading in field intensity, 
we make use of the maxima and minima of the curve E},(t), the effects of 
rapid and slow fading as well as interference and absorption fading, etc., 
and also certain statistical parameters. It should be noted that statistica 
studies of ionospheric radio wave propagation as the basis for compiling 
radio forecasts for the requirements of the technical services in commu- 
nications was confirmed only ten years ago. Among the statistical param- 
eters, we shall denote the average value, with respect to time, of E (the 
fixed component of the process) as practically equal, relative to the ave- 
rage over the total period, to mathematical expectations, and also to the 
so-called median value Eo, that is, the value of field intensity which ex- 
ceeds the instantaneous values during 50% of the time of observation. Fre- 
quently, the so-called quasi-maximum value is also used. This is the 
value of field intensity which exceeds the instantaneous values during 5% 
of the period of observation (some authors take 10% [3]). Through numerous 
observations it is an established fact that the average value is 35% of the 
quasi-maximum value [1, 2]. Other authors, by the way, point out that a 
small quantitative difference exists between the median and average values 
with respect to time [4]. Actually, with a Rayleigh distribution of the in- 
stantaneous values of field intensity, it is easy to establish the relation- 
ship between the average and median values as E = 1-063 Eo (i.e. the 
difference between the two values does not exceed 6.3%), while with a 
Gaussian distribution E = Ep [5]. For the so-called m - distribution, 
with certain exceptions, we also have either the exact equality E = Eo, or 
an approximate equality between these values (it goes without saying that 
our discussion concerns only radio wave propagation beyond the horizon) 
[6]. This shows that in a very large number of cases of practical, long- 
distance radio wave communication, the averaging operation gives not 
only the average value but also, to a first approximation, the median value 
of the electromagnetic field intensity. 

The literature makes reference to instrument-counters which auto- 
matically determine the various statistical parameters of a function of 
form E;,(t), since the direct determination of these parameters by means 
of operiocranis is very difficult [5]. These instruments record the in- 
stantaneous values of intensity for a specified interval of time. They in- 
dicate that portion of the interval in which the instantaneous values of in- 
tensity exceed a specified level or are confined to known limits. By means 
of the distribution curves obtained in this manner, the required statistical 
parameters are read off and calculated. These instruments are relatively 


‘The function of probability density of this distribution is determined in 
the following manner [6]: 
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where I'(m) — is the gamma function, the parameter m > 1[2 and Q = E2, 
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complex in their design. One version of this type of instrument is given in 
[7]. An electromechanical counter, connected in the plate circuit of an 
electron tube, starts to operate when a specified positive potential is ap- 
plied to the control grid of the tube. Ten such stages (10 tubes) are pro- 
vided; the value of voltage at which each of these stages begins to function 
is determined by the voltage which develops across its respective cathode 
resistor. The counter stages are cut off periodically by means of multi- 
vibrators, which makes it possible to make computations for a unit of 
time. 

Let us now proceed to certain theoretical observations. From Figure 1 
it is apparent that E),(t) is a continuous, random function of time which 
describes the fluctuations in field intensity. 2 It has been established, 
through numerous statistical data compiled during research on ionospheric 
radio wave propagation at various wavebands, that the instantaneous values 
of field intensity obey definite laws, to which we have already referred in 
references [2, 4, 5, 6]. It is reasonable to assume that these laws do not 
depend, on the one hand, on the season of the observations, and on the 
other hand, they maintain their strength independent of the duration and 
starting point of the observations.” This means that the function E).(t) is 
a stationary random function. As is known, a random process becomes a 
stationary [steady-state] one when all conditions creating it do not change 
in the course of time. Strictly speaking, one can hardly assume the ex- 
istence of such a process in the course of an unlimited interval of time. 
However, in practice it is quite satisfactory if these conditions remain 
constant during the total observation time and the time of the additional in- 
terval required for attenuation of the transient process [7]. 

In stationary processes, on the basis of the Erg hypothesis, the ave- 
rage value over the full range is identical to the average value with respect 
tOoLIMere e). 


+o be 
= 5 1 
E = | xp(x)dx=lim 7 J Feo ae 


—2o 


(1) 


Here, E is the average value of stationary function E},(t), p(x) is the 
function of the probability density of a random, stationary process, T — 
is the duration of the period of observation. 

It can be seen that the determination of the average value of field in- 
tensity may be obtained by calculating the integral in the second equality of 
(1). The proposed instrument is intended for the immediate determination 
of this value. From this it is possible, of course, to also determine the 
median value, providing we know (even if approximate'y) the corresponding 
distribution of the instantaneous values of field intensity. The literature 
gives valuable, practical instructions on this problem. For example, data 


2 With the condition that the source of the electromagnetic radiation and the 
receiving apparatus operate during the period of observation, and also with 
the condition that sudden absorption in the ionosphere does not occur. 


30f course, in practice this period is not limited. Thus, for example, a 
Rayleigh distribution in the meter and decimeter band is observed for a 
period of not more than 5 — 20 min. [5]. 
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is available concerning the distributions encountered in radio wavebands 
during a 24-hour period of observation, the distance between the receiving 
point and the radiation source, etc. [2, 3, 4, 5, 6]. Some of this data is 
contradictory [8]; in general, a complete investigation is far from realiza- 
tion. 

Taking into account the nature of the random process which is of in- 
terest to us, a simple electronic instrument was designed, the main com- 
ponent parts of which are a two-stage de amplifier with a relatively high 
input resistance and a de electric meter as the integrating element (Fig. 2). 
In its present form, our instrument is adapted 
for connection in the AGC (automatic gain con- 
trol) of a superheterodyne receiver with long, 
medium and short wavebands. Naturally, it 
may also be adapted to other wavelengths. 

The principal schematic diagram of the 
instrument is shown in Figure 3.4 The circuit 
for the de amplifier (with a common power 
supply) is similar to circuits encountered in 
so-called computing receivers. Under condi- 
tions of ionospheric propagation, the negative 
voltage at the output of the second detector (A) 
will vary from zero (minimum amplitude of the 
incoming signal to the receiver) to some nega- 
tive value (for maximum value of the signal). 
Thus, the changes in the resulting field inten- 
sity picked up by the receiving antenna will be expressed in the instrument 
through proportional changes in plate current of the second stage of the 
amplifier (from a certain minimum to a certain maximum value). At the 
same time, the rotor of the meter, which is connected in series in the 
plate circuit of this same stage, will rotate in step with the changes in cur- 
rent through its winding. At the end of the observation interval the number 


Figure 3 


*In connecting the instrument to the receiver, it is necessary to disconnect 


the lead through which the AGC voltage is applied to the grids of the first 
tubes of the receiver. 
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of rotations indicated on the meter will be a measure of the average value 
of electric current passing through it, i.e. a measure of the average value 
of field intensity at the receiving point. With potentiometer P9, the ampli- 
fication of the signal applied to the second tube may be varied in order to 
select the most suitable range for the plate current Ie]. Potentiometer P3 
is used to establish the initial operating point of the second tube (when a 
signal to the control grid of the first tube has not been applied). In order 
to decrease drift in the first tube, the well-known Miller circuit is used, 
which partially compensates for changes in heater voltage. Optimum com- 
pensation is achieved by means of potentiometer P,. The de voltage supply 
is provided by a conventional full-wave rectifier (the plate voltage of the 
right half of the first tube is regulated by means of a gas tube). The instru- 
ment is connected to a 220 volt power line through a ferroresonant regu- 
lator. To achieve greater stability in relation to drift, the instrument was 
provided with tubes which had already been in use for some time, and re- 
sistors with nominal power ratings twice as large, on the average, as 
called for in the design. As a result of the stabilization measures referred 
to above, a satisfactory reduction in drift was achieved — Fig. 4 (initial 
current through the coil of the electric meter Ig] = 30 ma; t = 0 at the in- 
stant of connection). From the figure it is evident that a preliminary warm- 
up of the instrument of not less than 1-1/2 hours is desirable. 


t (hours) 


Figure 4 


Since the maximum plate current of the second stage is in the order 
of several tens of milliamperes, whereas the nominal value of current 
through the electric meter usually amounts to several amperes, then it is 
expedient to disconnect the shunt mounted in the meter, with the result 
that the plate current of the second tube will now pass directly through the 
rotor winding of the meter. Knowing the electrical data of the meter, it is 
easy to derive the following equation for the current passing through it: 


60 Ry + lan MOS 


1000 Rsh ie IC 2 («= 1000 Rsh ) 
60 Ry»tr U T Fu 


do 


(2) 


Here I) (ma) — is the current through the rotor, Rsh (ohms) — is the 
de electrical resistance of the meter shunt, r (ohms) — is the de electrical 
resistance of the rotor winding, k' (watt-sec/rev) — is the meter constant, 
U (volts) — is the nominal meter voltage, n (rev) — revolutions of the rotor 
as recorded on the meter, T (min) — period of observation. : 

An important characteristic of each meter is its error curve, Figure 5 
shows its approximate form. As can be seen, the error is practically con- 
stant beginning with a known value of current (this value is about 30 ma in 


the figure), and it may 
[* readily be taken into account 
in calculating the number of 
revolutions. 

Figures 6 and7 show two 
curves of I,] = @(-Uin) for 
the dec amplifier: for low and 
mda) high currents. A feature of 

this amplifier is that it op- 

erates with a certain initial 
plate current (second stage). 

Thus, the following advan- 

tages result: 1) a decrease 

in the effect of the curvature 
: at the lower part of the plate 
Figure 5 current — grid bias charac- 
teristic of the second tube, 
and 2) compensation of the inherent error in the meter. The error of the 
instrument does not exceed 2-3% on the average. It is desirable to deter- 
mine the number of revolutions to an accuracy of a few revolutions, and 
T — to an accuracy of seconds [Formula (2) can be modified readily], 
which practically eliminates errors in the determination of these param- 
eters. 

If the number of revolutions and period of observation are known, then 
from formula (2) it is easy to calculate the average current I,). For con- 
venience, a family of curves of I,] = © 1(n) (parameter T) may be plotted; 
in this case I,; may be computed directly from the corresponding T - curve. 

The receiver to which 
the instrument is connected 
is calibrated beforehand: the 
relation E = @9(-U;,) is 
known (E — is the field in- 
tensity at the receiving 
point). Therefore, the un- 
known average (and corre- 
spondingly, the median) val- 
ue of the observation may be 
found in the following manner 
(Figure 8): calculate or read 
off the current I,] =Io, then 
mark this value (Ig) on the 
ordinate of Ig) and finally, 
follow the procedure indica- 
ted by the arrows in Figure 8. 
As a result, the unknown val- 
ue Eg may be read from the 
second ordinate (E). 

In testing the instrument 
a series of observations were 
carried out at a frequency of 
155 ke — Brashov radio 
transmitter (distance from 
receiving point — 380 km). 
For this purpose, a re- 
Figure 6 corder was connected to the 
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instrument. From the re- 
corded data and usual methods 
of calculation, the correspond- 
ing median values were deter- 
mined. The latter were com- 
pared with the average values 
indicated by the instrument. 
During the initial determina- 
tion of median values, inac- 
uracies may be permitted 
since the recordings are 
frequently quite complex. Be- 
sides, it is difficult to reach 
exact conclusions regarding 
the law of distribution of the 
instantaneous values of field 
intensity for a limited number 
of observations. Of course, 
this was not one of the pur- 
poses of the trial tests. Nev- 
ertheless, on the basis of the 
values in Table 1 for relative 
difference between the ave- 
rage (E) and median values 
(Eg), it may be stated, gen- 
erally, that the observed 
distribution of the instantan- 


eous values of field intensity at this wavelength is of the Gaussian type. 
During observation No. 9, this distribution seems to approach the Ray- 
leigh type, and for observations Nos. 8, 10 and 11, it is close to an m-type 


(mV/t) 


A (ma) 


Figure 8 


distribution. The inaccuracies 
in observations No. 6 and 12 
can be explained by insuffi- 
cient preliminary warmup of 
the instrument (these were 
initial observations). The 
tests, which were carried out 
on September 7, 1958, show 
that for different periods of 
observation, different dis- 
tributions may be obtained. 
The instrument will op- 
erate satisfactorily up toa 
fading frequency of about 
2 cps. As was reported at 
the MKK conference of 1953 
in London, the fading fre- 


quency may reach 10 ke. However, this will not preclude the wide applica - 
tion of the described instrument since very rapid fading of this type is, in 
practice, of less significance in conventional statistical studies. Also, they 
are not reproduced by the recording inks of recording instruments which 
are commonly used for surveys on radio wave propagation. 

The instrument has good operating qualities. Maintenance is also very 
simple. It is recommended that the recording apparatus be operated simul- 
taneously with the instrument — to monitor the field, particularly in case 
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Table 1 


Fo (mV/m) 


4,1X 1958 


20152115 
21152145 
21452217 
21 1a NF 


9/1X 1958 


10/IX 1958 


11/1X 1958 18°0— 18% 
1830—1900 
19302300 


of ionospheric disturbances, and also for extracting additional useful in- 
formation. A shortcoming of the instrument is that the determination of 
the median value (from the average) becomes impossible when the desired 
signal is lower than the noise level for a period of even less than 50% of 
the total period of observation, despite the fact that this value does exist 
in this case. 
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SYNTHESIS OF FILTERS 
WITH LINEAR PHASE CHARACTERISTICS 


A.F. Beletskiy 


Several methods for solving the problem of approximating the transfer 
functions of filters with linear characteristics as well as, in the general 
case, with given phase (shift) characteristics, are examined. It is assumed 
that the synthesized filter will be a reactive four-terminal network. 


In recent years it has become increasingly necessary in communica - 
tion engineering to use filters which, in a specific frequency band, have 
phase shift characteristics approaching linearity in one sense or another. ! 

The current article is devoted to the examination of several of a num- 
ber of possible methods for synthesizing electric filters with linear phase 
characteristics by means of their operating parameters. It is assumed 
that the synthesized filter will take the form of a reactive four-terminal 
network. 

The problem of synthesizing a filter of one of the various types to be 
discussed below shall be considered solved if the function expression the 
filter transfer function satisfies the conditions of the problem, as well as 
cunditions for its physical realization. This implies that we shall examine 
all practicable methods for approximating the transfer function of filters 
with linear phase characteristics. In order to determine the appropriate 
circuits and values for the elements comprising the possible variants of 
these filters, we shall make use of the known methods for the synthesis of 
reactive four-terminal networks by means of a given transfer function. 


lHenceforth, for the sake of simplicity, filters of this type will be referred 
to as filters with linear phase characteristics. 


The results obtained below may also be broadened to include the syn- 
thesis of filters with nonlinear phase characteristics. 


PHASE CHARACTERISTICS OF REACTIVE FOUR-TERMINAL NETWORKS 
When a reactive four-terminal network is connected between purely 
resistive loads (Fig. 1), its effective propagation constant g, effective at- 
tenuation b, effective phase a and effective transfer function S(i@)/are re= 
lated by the following relationship 
g=b -ia=InS(io), (1) 


where 


S (io)= a VA __ 4 (iw)? + ay (lo). an ; 


9 Re (es) 7 ee (1) eh 


Let us recall [1, 2] the properties of the function 
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The Hurwitz polynomial v(p) 

appears in its numerator, i.e. a 

polynomial, with real coefficients, 
AL 5 all roots of which are found to be in 
. RI tu, the left hand p plane. 
E ap The denominator of the function 

is, in the general case, an even or 
odd polynomial f(p), all coefficients 
of which are also real. It should be 
noted that, in particular cases, the 
polynomials v(p) and f(p) in Eq. (2) 
may have a common multiplier ~(p). The latter, on the basis of the prop- 
erties of the polynomial v(p), is the Hurwitz polynomial. Under these con- 
ditions the polynomial f(p) will be neither an odd nor an even one. However, 
it may always be reduced to one of these, without disturbing the properties 
of the function S(p), by multiplying its numerator and denominator by the 
Hurwitz polynomial u (p). 


Finally, with a change in p along the imaginary axis, i.e. with 
p= iw 


Figure 1 


|S (iw)| > 1. (3) 


The effective phase of a reactive four-terminal network is equal to 


a = arg [S (iw)] = arg Freee 


(4) 
Since f(i«) can only be an even or odd polynomial, then it follows from 
the latter relationship that the effective phase of a four-terminal network 
will be determined — to an accuracy of the component r7, if f(iu) is an 
even polynomial, or of the component (r — 0.5)7, if f(iw) is an odd 
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polynomial — by the argument of the Hurwitz polynomial v(iw), i.e. by the 
argument of the numerator of the polynomial S(iw) [3]. Actually, as w 
changes, the polynomial f(i#) — which has either purely real values when 
f(iw) is an even polynomial, or purely imaginary values when f(iw) is an 
odd polynomial — may change in sign when real roots in odd multiples are 
present. A change in the sign of the polynomial also leads to a jump in 
phase by an angle of 7 radians. It should be noted that for each jump in 
phase there is a corresponding pole in the transfer function S(iw) and a 
discontinuity in attenuation of the filter. 

In accordance with the above statements 


1 eters? Gate) mae mN 


T 


a=arc tg — 


fl wie) z(t) = © 9) (5) 


It is important to note that the poles of the function S(p) for a four- 
terminal network of the minimum-phase type may be located only on the ima- 
ginary axis ofp, i.e. whenp=iw, including the values w=Oandw=%, This is 
a direct result of the even (odd) polynomial f(p) of the denominator in 
Eq. (2). Therefore, in reactive four-terminal networks of the minimum- 
phase type, if they are not polynomials or semi-polynomials [2], then 
jumps in phase are unavoidable. 

In conclusion, it should be noted that if we substitute the inequality 
[S(Gi@)| > 0 for the inequality (3), then all of the above statements hold 
true not only for the effective phase of a reactive four-terminal network 
but also for phase shift in current, voltage and electromotive force. 


POSSIBILITIES OF PRODUCING FILTERS WITH 
GIVEN PHASE AND ATTENUATION CHARACTERISTICS 


In the general case, the number of variable parameters for the trans- 
fer function of a reactive four-terminal network, which is a rational func- 
tion of type (2), does not exceed 1.5 n + 1 for even values of n, and 
1.5 n + 0.5 for odd values of n. 

For approximating a given relationship between filter phase and fre- 

- quency, no more than n parameters should be used since the filter phase 
is determined by the polynomial v(iw) of function (2). The given relation- 
ship between filter attenuation and frequency may be approximated by 
using all the variable parameters of this same function. In the latter case, 
the filter phase characteristic is not controlled, which is typical for stand- 
ard methods of synthesizing filters on the basis of their operating param- 
eters. In the general case for approximating a given relationship between 
filter attenuation and frequency, the number of variable parameters used 
may be equal to the difference between their total number and the number 
’ of parameters used for approximating the given relationship between filter 
phase and frequency. 

If the conditions of the problem specify that phase distortions in the 
filter within the boundaries of its pass band must not exceed a certain 
permissible value, then these requirements may always be satisfied by 
using a Hurwitz polynomial v(p) of as low a degree as desired, beginning 
with n = 1, while the requirements for the relationship between attenua- 
tion and frequency may be fulfilled, usually, only by functions of the form 


*tg = tangent (Transl. ) 
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(2) of relatively high order. Consequently, a problem arises in regard to 
the determination of the most suitable relationship, between the number of 
variable parameters used for approximating the given relationships of filter 
phase and attenuation, which will result in a minimum number of elements 
in the given or selected filter structure. A similar problem arises in the 
general case when the given phase relationship of the filter within the bound- 
aries of its passband is nonlinear. 

In the above-indicated general treatment, the problem of systhesizing 
filters on the basis of the given requirements for their phase and attenua- 
tion characteristics is seldom solved analytically, with the exception, of 
course, of its simplest variants. 

However, it is always possible to select a Hurwitz polynomial v(p) of 
sufficiently high degree so that the requirements in the relationship between 
phase and frequency of the filter in its passband will be satisfied. 

By further selecting a corresponding polynomial f(p), the degree of 
which must not exceed the degree of the polynomial v(p), itis then possible 
to fulfill the requirements for the relationship between filter attenuation and 
frequency. At the same time, the filter phase characteristic, which was 
either found or selected in the initial stages, will remain unchanged. 

The solution of the relationship between filter attenuation and fre- 
quency is a function of frequency of the form 


_ 1) v(lwyo(—te) Vo (iw) vo (— te) 
lig up ayy GR & te DeamiG (reels (6) 
The resulting solution will not be an optimal one in the general case 
for reasons indicated above. A numerical determination of the loss in the 
number of filter elements is not possible, mostly because of the absence 
of an optimal solution of the problem. However, this loss will seldom be 
a significant factor if sufficiently rigid requirements are set for accuracy 
in reproducing the given relationship between phase and frequency. 
The given phase-frequency characteristic may be approximated by 
again using a Hurwitz polynomial of the form 


U(p) = vf (p). 


Corresponding to this polynomial, we have 


(7) 


Sime = [ye Us het jw)]? [ wg(iw) 9 
J (iw) J (iw) Vo (— iw) | (8) 


The product 
Vo (w?) = ty (i@) v (— io) 


is, obviously, an even polynomial with real coefficients, where 
V, (w?) > 0 


for any real values of w. 


A transfer function of the form (8) will have a corresponding phase 
characteristic 


Uo (iw) =p (— iw) 
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which differs from the phase characteristic of a network composed of 


uniform phase circuits only by the component r 


T 
2 

The attenuation factor in this case may be determined from the func- 
tion 


mai (w?) 


== iba —, 
If (¢@)| (10) 


Here, f(iw) represents the modulus of an even or odd polynomial, all 
coefficients of which are real. 


LOW FREQUENCY FILTERS WITH LINEAR PHASE CHARACTERISTICS 
AND WITH ISOEXTREME ATTENUATION CHARACTERISTICS 
IN THE ATTENUATION BAND 


Let us solve the problem of synthesizing a low frequency filter in 
which: 

a) the transfer function is expressed by a function of form (2); 

b) in the frequency band 0 < W< Wog, either the deviation of the 
phase-frequency relationship from linearity, or the deviation in delay 
time, from a certain fixed value, does not exceed a value which is con- 
sidered permissible; 

c) the attenuation for u!= 0 is equal to zero; 

d) in the frequency band «> Wop, the smallest attenuation values 
will occur when the remaining equal conditions are the maximum possible. 

Proceeding to a solution of the problem, let us select a Hurwitz poly- 
nomial v(p) of degree n = 2! , the phase characteristic of which satisfies 
condition (b) of the problem. 

As a polynomial of this type we may use, for example, a Bessel 
polynomial [4, 5, 6] of a corresponding degree. 

Then, in accordance with (6), the attenuation of the synthesized filter 
will be determined by the function 


2 y2 
(2) eee ee en bo aie ep del 


Vo (io) (10) VV cae ae Byes feo ae (11) 


In this expression, the degree of the polynomial 


Ff (iw) = Boo” + Bw? +...- 1 (12) 


is taken as high as possible, which corresponds to the maximum possible 
number in the selected polynomial 


V (w2) = 0,0" + a0 +...4 1 te 


2Filters with linear phase characteristics which correspond to the par- 


ticular case when q = 1 are examined in reference [3]. 
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of variable parameters which are used in approximating the attentuation 
characteristic of the filter. In addition, expression (11) is standardized 
so that condition (c) of the problem under examination will be satisfied. 
It is also necessary that in the interval & > Wop, the maximum de- 
viation of function (11) from zero will be as small as possible. 
Using methods employed in the theory of highest approximation of a 
function [2], we find that 


: oP: 
a? — Oa ba ( “3b a e —a?)? + 4p o, 
cos y) arc COS 
F (w?) = r=1 


V (eet —eyt ee a 
I 


2 2 2 2\2 222 

— &,— wo, +P, — 2 + 4prc 

cos ) arc COS ae Vv ( * i a = 
Pr 1S, 


(14) 


r=1 


Here (-py + idy) are the roots of the polynomial v(p) with a positive 
imaginary component. 

All null points of function (14) lie inside the interval of its least de- 
viation from zero wo, < w < ©. Consequently, all discontinuities in the 
filter attenuation 


b= —In|F (@?)| (15) 


lie in the attenuation band of the filter, which condition is satisfied by a 
four-terminal network of a minimum-phase type. The total number of dis- 
continuities in filter attenuation will be equal to the degree of the poly- 
nomial f(iw) with respect to the variable w2, i.e. equal to /. 

From the latter relationships it follows that within the limits of the 
filter attenuation band 


: 2 V 2 2 .2)2 2 2 
b+ ( ob + P- —<;) + 40, 9, 
b>by=ch } Arch Rae 5 (16) 
r=1 
Here, the value of the guaranteed filter attenuation in its attenuation 
band bg corresponds to the sign of equality. This lowest value of filter at- 
tenuation passes at/ + 1 frequencies, including the frequencies Ww = Wob 
andW=°. 

In the frequency band0< w< Hop: the filter attenuation will increase 
monotonously from the value b = 0 to b = by, providing the polynomial v(p) 
is a Bessel polynomial. A similar type of change in filter attenuation within 
the limits of the indicated band is also characteristic for a number of other 
types of Hurwitz polynomials. 

Figure 2 shows the dependence of filter attenuation and delay time on 
frequency, which was calculated by the method presented above for the 
particular case when n = 6 and wop = 5. In the example, a Bessel poly- 
nomial of sixth degree is used for the polynomial v(p). 

For each given type of Hurwitz polynomial, it is possible to calculate 
the relationship between the guaranteed filter attenuation bg and the de- 
gree n of the polynomial, and the value of the limiting frequency of the 
quaranteed attenuation band woh. This relationship, which was calculated 
by A.T. Lebedev using Bessel polynomials, is shown in Figure 3. The 


*ch = cosh (Transl. ) 
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Figure 3 


phase characteristics of filters in this case do not differ, with respect to 
the above considerations, from those expressed in Bessel polynomials. 
The corresponding relationship between delay time and frequency is shown 
in the same figure. 


GRAPHICAL METHOD FOR CONSTRUCTING THE 
ATTENUATION CHARACTERISTICS OF FILTERS WITH 
GIVEN PHASE CHARACTERISTICS 


Function (6) may be transformed into the form 
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If w,, are real numbers, then the relationship between frequency and 
any one of the components within the summation mark of the latter ex- 
pression may be constructed on a logarithmic frequency scale by means 
of parallel transformation, using the same pattern in each case. The 
shape of the pattern is determined from the relationship 


6, = — (u+In2sh|u)), (18) 
where 
u=In| -—— 


The pattern is shown in Figure 4. the first component of (17) is de- 
termined, to an accuracy of the fixed component, by the polynomial v(iw) 
which is initially selected in accordance with the requirements of the fil- 
ter phase characteristic. The dependence of this component on frequency 
may be shown to the same scale as in (18). 

- By using the described pattern and changing the number and location 
of the attenuation discontinuity of the filter, it is then possible to fulfill 
the requirements of its attenuation characteristic in the stop band. In the 
pass band the attenuation characteristic for low-frequency filters will be 
similar to that shown in Figure 2. 

Naturally, the described pattern is suitable for the construction of 
attenuation characteristics of band-pass filters as well as of low-frequency 
filters. In order to increase the attenuation of such filters in the low fre- 
quency region, it is permissible to insert the attenuation discontinuities 
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3 


Figure 4 
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for (ea = 0. On a logarithmic frequency scale, for each such attenuation 
discontinuity of the filter there will be a corresponding straight line with 
a slope of 0.69 nepers per octave. 


6 nepers 


Figure 5 


Amplitude distortions within the limits of the pass band of a low- 
frequency filter may be reduced if we permit the existence of imaginary 
roots in the plynomial f(iw). In this case, it is necessary to use a pattern 
in which 


6,= —(u+1n2chu). (19) 
This pattern is shown in Figure 5. 


CORRECTION OF FREQUENCY ATTENUATION CHARACTERISTICS OF 
FILTERS WITH LINEAR PHASE CHARACTERISTICS 


In a number of cases, the attenuation characteristics for filters of 
the type examined above prove to be unsuitable because of significant am- 
plitude distortions in the pass band. 

In such cases, some of the variable parameters in the polynomial 
f(iw) of a function of the type given in (6) must be used to control the at- 
tenuation characteristics in the pass band. The latter causes, as a rule, 
the appearance in the plynomial f(iw) of roots outside of the real axis of the 
variable w, and correspondingly, a transition to structures of the non- 
minimal phase type. 

As yet, a rigorous solution of the problem, by methods employing the 
- theory of highest approximation of the functions, has not been attained for 

the case of filters with linear phase characteristics and with isoextreme 
attenuation characteristics both in the stop band and in the pass band. 

The attempt to use, in this case, solutions which are not optimal ina 
Chebyshev sense will, as calculations show, result in filters with an un- 
realistically high number of elements. 

The simplest of all, for our purposes, is to equalize the attenuation 
characteristics of the filter in its pass band with the aid of these same 
functions of form (6). The values of their variable parameters may be de- 
termined most conveniently by the interpolation method. Then the trans- 
fer function of the filter with an equalized attenuation characteristic will 
be equal to the product of the transmission constants of the equalizer 


. 


- 
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network and of the original filter, which has an unequalized attenuation- 


frequency relationship in the pass band. : 
As an example, the solid lines in Figure 6 represent the attenuation 


0 } sh i lw 
U Zen Oh 45. e268 B 


Figure 6 


characteristic and delay time of a filter calculated by the above methods 
for ‘the particular case when N = 8 and Wop, = 7. On the same graph, the 
dotted line represents the attenuation characteristic of the filter after 
its correction by means of the function 


|_ (3) 


6, = In »————_ 
i lagu’ + ayo? + ay\’ 


where vi) is a Bessel polynomial of seventh degree. The latter intro- 


duces practically no phase distortions in the given frequency band 0 q¢w <5. 2. 
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STRUCTURAL STABILITY 
OF A LINEAR FOUR-TERMINAL NETWORK 


V.Z. Nayderov and Yu. L. Simonov 


The basis for a method of increasing the stability of a linear four- 
terminal network by correcting its internal feedback is presented. A struc- 
ture for the four-terminal network which is stable for any load on its ter- 
minals is shown to be attainable. 


One of the causes for operating instability in linear amplifiers is in- 
ternal feedback in the amplifying element of the circuit (electron tube, 
transistor, etc.). In designing amplifiers, this feedback must be taken 
into account, or else its effect must be compensated for by additional ex- 
ternal feedback (various neutralizing circuits). 

The proof of one of the theorems on the stability of an active four- 
terminal network is given below. 

Using a transistorized tuned amplifier as an example, a method for 
increasing the stability of linear amplifiers, based on the use of the proven 
theorem, is examined. 


STABILITY OF A LINEAR FOUR-TERMINAL NETWORK 


Theorem. A non-autonomous, active, linear four-terminal in which 
the transfer function is determined by the expression 


OQ Yu+¥n' (1) 


1Regarding the problem of stability of linear four-terminal networks, the 
editors propose to review the results of discussion on the article by A.A. 


Kulikovskiy (Elektrosvyaz' No. 12, 1959). 


2This condition signifies that the phase shift between Ug and Uy equals 
180°, if we ignore the additional phase shift introduced by the passive 
elements of the four-terminal network. 
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I 1 will be absolutely stable for any 

t passive loads on its input and out- 
put terminals (Yo and Yy,) provi- 
ing it is stable when the input and 
output terminals are shortcircuited, 
and if the product of the character- 
istic conductances Y,,Y9; are real 
numbers. 

The physical processes in a non-autonomous linear four-terminal 
network (Figure 1) which satisfies condition (1) and which is loaded at the 
input and output terminals by the conductances Yj and Yj, respectively, may 
be described by the following system of equations, written in matrix form 


a 3a Yio Wedel: | 
= Ney — (Yo + Y) U, !, 


Here, Y4,...Y2.are the characteristic conductances of the four-ter- 
minal network. 

In reference [1], it was shown that a four-terminal network will be 
absolutely stable if its amplitude-phase characteristic 


Yi2 Yor 


8K = =a+ib 2 
: (Yut Yon + Yd (2) 
does not encompass the point BK =1. 
Assuming: 
Yun ze Yo= Bu +e id, — Vi = £12 + 1 dj2, ¥3; = £ + i Bo, 


Yoo + Yi = B22 + i bos, 


considering that 


Yy2 Yai = Re (Y19 Yo) = —(812 821 — P12bg1), and sepa- 
rating the real and imaginary parts in (2), we obtain the following expres- 
sion for a and b: 


ae IY ol? (611402 — B11 B22) 
812 IY Yo9i? ; 
he Se ial? (611829 + 609 £11) 
812 |" 1Y 2012 : 


(3) 
The requirement for stability is met when a < 1 for b = 0; accordingly, 
b,=— bo» fu, 
&n2 (4) 
Substituting (4) in (3), we obtain 


y 


2 
Yn 


oS = 811 821 
812 £22 


(5) 


wk The right part of this equation is negative for any real values inserted 
in its terms; consequently, a < 0. The BK diagram excludes not only the 
point BK = 1, but also the point RK = 0. The four-terminal network will be 
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stable with a sufficiently large stability margin for any values of Yo and 
YT: 

In the particular case of a transistorized tuned amplifier, the proof of 
the theorem may be carried out on the basis of the stability conditions ob- 
tained in reference [4], 


Roe (Xe Wi2¥a 
(822 Roe + oe) [1¥2¥o1! + Re (¥12¥o1)] 


2g1Roe < 2n?, (1 —k,). 


Assuming in this inequality, in accordance with the theorem, that 
Xe (Y12Y9;) = 0, it is easy to see that the tuned amplifier will be stable for 
any values of njy, Noyt and Rog. 

The proven characteristic may be used to increase the stability of 
linear amplifiers, for example, transistorized tuned amplifiers. 


INCREASING THE STABILITY OF A TRANSISTORIZED SELECTIVE 
AMPLIFIER WITH A COMMON-EMITTER CIRCUIT CONNECTION 


The characteristic conductances of a transistor in a common-emitter 
circuit are expressed most frequently in the following form [2]: 


ee Alt i ORS 
ri 
1 F 1 . 
a —— +ioC,.= — (1 +i), 
rio ry 
ve 1 sl 1 
try + iw Loy ra (1 + i %9)) ’ (6) 
Yo9 = +10Cy., 
T2 
where 
iE 
19=M2C yo, t= re . 

y It is not difficult to show 
that, in this case, the condi- 
tion Y,;9¥Y21 = Re (Yy2Ya1) 
takes the form 

So: “ore (7) 
Figure 2 ~ Jn order to achieve this, 


an additional passive conduc- 


tance Y is connected to the transistor (Figure 2). 
The matrix for this combination will be 


Fire all Yi ey aa 
— Vine —Y ne, == Ne an — A 


Ly= | i 


The matrix elements in (8) may be written in a form analogous to (6): 


1 


re 


Yye= V+ Y= +i0Cire 


ol 


1 1 


—VYre=—VYat+ Y= + i@Cpre= (1 + 1 %2¢) 
: rye rie 
1 1 
Y.,e= Y,., — Y= —— ; 
ga = ryetiwLloe roe (1 + 1 tae) 
1 zZ 
Vee a Y= — + ioC,, 
22€ 29 + Fal 22e (9) 
where 
= 21e 
Tize=SiveC xe) t21e = —" 
rae 


For an equivalent four-terminal network, expression (7) may be re- 
written in the following form: 


“126=— “ele (10) 
Table 1 shows the relationships connecting the components of the char- 


acteristic conductances Yj9¢ and Yo;,, respectively, with the components 
of the conductances Yj9 and Yo, for various forms of the conductance Y. 


Table 1 


ry2(Cy2 + C) 
ra 


(1 + w2rayt21C)? + 


Tay? (roto Sa 
ao +( Lavon + wr C2 


ta +ra1C (1+ w2r5,) 


By substituting the values T 95, and T 21e from the table into (10), it is 
easy to obtain the following expressions which connect R, L and C with the 
parameters of the triode frequency: 


r 
aie 7h + 023) 
1 


Me (for % 4. >), 


1—4 (1 + «28 
r 


L= 2 » (for 1. >t), 
02(C + Cc; = wa ) 
ry2 
ae (11) 
Cc = 7 a (for So4 = X12) 
VS eae + 0223) 
where C, — is the self-capacitance of the inductance L. 


Calculations show that in the frequency band where it is usually nec- 
essary to pay particular attention to increasing the stability of selective 
amplifiers, the following inequalities must be adequately fulfilled: 


C72 (Ip ot)< 1, “(1 +0?) <1. (12) 

then [ee ee d , C=C, (= =p 
ah Selle A 721 TH2 13 
12 12 2 (Cx -++ Cy. re ) ( ) 


As is obvious from (8) and (9), when a resistance R is used as the 
element for increasing the stability, the triode output resistance Rout may 
decrease considerably because of the increase in Yo9¢ in comparison with 
Y92. This leads to an undesirable shunting of the oscillatory loop in the col- 
lector circuit. Therefore, an inductance should be used in place of the re- 
sistance. 

However, from constructive considerations for Tj9 > 7% 1, it is more 
expedient to use capacitance instead of inductance as the element for in- 
creasing the stability. To facilitate this substitution, the amplifier circuit 
should include a phase-inverting autotransformer or transformer(Fig. 3). 

It is not difficult to show 
that the equivalent characteris- 
tic conductances of this circuit 
may be determined from the fol- 
lowing relationships: 


Vi2e= Vin + io @ C, 
Pi 


ape IE 
Y,,.=Y sea 
21e a + 10 P, C, 


where p, and po — are the coefficients for the oscillatory circuits con- 
nected into the collector and base circuits, respectively. 

Determining, from this, the expressions for 712eand791e, as was 
done in (9), and substituting them in (6), we find the following formula for 


the capacitance C: 
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This formula may be simplified if we disregard the small terms in 
accordance with (12). 
We have 


C=C), sail = el 


erp? 712 (14) 
The circuit of Fig. 3 is analogous to an amplifier circuit with neu- 
tralization of the y - type, but the value of the capacitance, determined 
from (14), is lower than that calculated by formula (11) in reference [3] 
by the magnitude 
Crete ale 
P2442 


Thus, the circuit (Fig. 3) may be used successfully in place of a neu- 
tralizing circuit of the y - type in which the neutralizing network is com- 
posed of resistors and condensers. A substitution of this type is feasible 
only with weak coupling between the triode and the input and output oscil- 
latory circuits, or with relatively low Q - factors for these circuits, since 
such shortcomings as interaction between the tuning circuits have not been 
completely eliminated. 


EXPERIMENTAL RESULTS 


The experimental verification of the theorem which was proved above 
was carried out on a tuned amplifier (Figure 4). 


A-ChM-2 A-ChM-2 
amplifier (K = 50) SI-1 
q a) 


Figure 4 


The circuit makes it possible to carry out measurements of: a) the 


U2 
U , b) the stability 


1 4 
factor ky ~ ica (Here U, — is the resonant voltage across the LyCy 


resonance amplification factor by the voltage Ko = 


circuit with connection Py closed; UL — is the maximum resonant voltage 


across the same circuit with connection P, open and a corresponding 
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detuning of the Lg, Co circuit), c) the bandwidth for each of the oscillatory 
circuits. 

The following types of transistors were used in the circuit: P-15 @re- 
quency range 100-300 ke), P-403 and P-405 (frequency range 1-5 Mc). 

In all of the measurements, the values for the connection coefficients 
p; and Py were selected so that the frequency bands passed by the oscilla- 
tory circuits would be approximately equal to each other. The measurement 
results for ky and Ko for the P-15 triode are shown in Figure 5). 


50 
100 40 180 220 = 260 = 300 m0 800 0 220 260 300 


Figure 5 


The numeration of the curves corresponds to the cases: 1) Regular 
amplifier (p; = 0.012, Po = 0.2). 2) Regular amplifier (py = 0.025, 

Pg = 0.38). 3) Amplifier with increased stability. An inductance L, cal- 
culated by formula (13) for a frequency of 150 kc, was connected across 
the terminals 1 - 1 (p, = 0.025, po = 0.38, L = 24.6 Mc [Sic]). 4) Ampli- 
fier with increased stability. A resistance R, calculated by formula (13), 
was connected across the terminals 1 - 1 (py = 0.038, po = 0.25, R = 36 
kohms). 

From an examination of the graphs, it may be concluded that the pro- 
posed method for increasing the stability permits a substantial gain in 
amplification of tuned amplifiers without decreasing the stability margin; 
in this connection, the use of an inductance L permits an increase in sta- 
bility in a narrow frequency band, whereas a resistor R serves the same 
purpose in a wide frequency band. Similar results were obtained with the 


triodes P-403 and P-405. 
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USING TRANSISTORS 
A QUARTZ MASTER OSCILLATOR 


S.L Kunina 


A master oscillator using a P-403 transistor and a quartz crystal be- 
tween collector and base is investigated in the 0.75 to 5 Mc frequency 
band. The effects of operating conditions and circuit parameters on fre- 
quency instability during changes in temperature and line voltage are studied. 
An engineering method is developed for carrying out calculations on the 
master oscillator. 


INTRODUCTION 


In references [1-3], which were devoted to an investigation of crystal 
oscillators with junction-type semiconductor triodes (PPT), the results of 
measurements of frequency instability were obtained for circuits with the 
crystal inserted between collector and base. Expressions were worked out 
for calculating amplitude and frequency of auto-oscillations at the limits of 
the self-excitation region. However, none of these references presents a 
detailed investigation of frequency instability under steady-state conditions, 
or recommendations in regard to the selection of operating conditions and 
circuit parameters. 

The purpose of our work was to investigate the properties of a transis- 
torized master oscillator, with a crystal connected between collector and 
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base (Fig. 1), at frequencies of f < 0.05 ie (where the inertial effects of 
the transistor on the operation of the oscillator are small), to study the ef- 
fects of operating conditions and circuit parameters on frequency insta- 
bility, and finally, to develop methods for engineering design which take 
frequency instability into account. 

For analogous oscillator circuits using electron tubes, similar prob- 
lems were solved by V.P. Zhu- 
khovitskay in references [4, 5]. 
In our work, the analytical meth- 
od developed in [4, 5] was used, 
taking into consideration the 
properties of the static charac- 
teristics and high-frequency 
parameters of the transistor. 

The design relationships 
were obtained in a system of 
Y-parameters for a common 
emitter circuit. 

These investigations were 
carried out with the use of a 
fragmented-line approximation 
for the static collector current 
characteristic, fragmented- 
parabolic approximation for the 
base current characteristic [6] 
and, in addition, the following 
expressions were used for the 
relationships between the admittances and parameters of the triode [7]: 


Figure 1 


Vi, = Sy) = So + iwCoe (1) 
Vn. = iwC,. (2) 
¥,,=S=Se™, (3) 


Here, Co, and Cj, are the input and output capacitances of the triode, 
So — the input conductance, S, ©, — the modulus and phase, respectively, 
_ of the collector current characteristic. 

We did not obtain an expression for the admittance Yy9> nor have we 
in Eq. (2) taken into account the active component of the output conductance 
Zke- This is explained by the fact that for f < 9.05 fg, under oscillating 
conditions (with a small load resistance), the effects of Yj and gye¢ on con- 
‘ ditions and frequency of self-oscillations are very small [7]. 

At frequencies of f < 0.05 fy, the phase @, of the characteristic is 
usually small and Eq. (3) may be replaced by the following approximation: 


Y= S~ S(1+4 ig,). (4) 


bie — is the frequency at which the modulus of the current transfer factor 


of the triode in a common-base configuration will drop V2 times. 
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The parameters of the triode C;,,, Coe, Sq; § and Os depend on fre- 


quency, collector voltages E,, current I, and temperature t°. 


At any frequency (whenf < 0.05 f,,) the expressions (1-4) may be used 
for analyzing the processes in master oscillators not only at the threshold 
of self-oscillations, i.e. when the fundamental voltage U, ~ 0, but also for 
investigating stationary conditions with high amplitudes for U and cutoff 
of collector current. The difference is only in the method used to determine 
the transistor parameters. In the first place, at any frequency for each of 
the values Ej,, I, and t’, the local parameters Cke, Coe, So. §, and ,may be 


measured for amplitudes U 0 +0. However, when investigating stationary 
conditions, quasi-linear (average) values of the parameters Chey: Coe1> 


So ,> 51 and may be calculated for each frequency and temperature in 


siZ 


relationship to Che: Coe: So: Ss, O. (E,., hh) by an approximate method pre- 


sented in reference [8]. 

The capacitances Cke1 and Cge1, which are connected in parallel with 
the capacitances C, and Co of the circuit (Fig. 1), are for convenience 
considered to be part of C, and Cy. Then, in analyzing stationary condi- 
tions of the oscillator, it is sufficient to take into account only the following 
triode conductances: 


Yai ey sy Sov (5) 
Yo, 1, =S, = S, (1 +i 9,)). (6) 


In accordance with [8], 5; and So; are also calculated as average 
plate and grid current transconductances of tube [9]. In using the approxi- 
mations, the characteristics Sy and S01 are determined from the expres- 
sions in [8]: 


Seo a) Cin) 
So, = aU, 1\?) (65). (8) 


Here, a — is the coefficient of the parabola approximating the base 
current characteristic; 6, 9 9 — are the angles of cutoff of the collector 
and base current (usually in semiconductor oscillators 6 = 86 [6]; ¥,(9), 
y,-*) (8) — are the resolution factors for a cosinusoidal pulse of the 
first harmonic with line-fragment and parabolic-fragment approximation 
of the characteristics. In reference [6] tables are given of the values of 
yf?) (8) and y (2) (6)2 for various 6 

1 0 ; 

The expression for calculating the phases of the quasi-linear trans- 
conductance (gs, was not worked out by us since research shows that for 
f < 0.05 f,, this factor has no influence on oscillator operation. The 
method for calculating the changes in @,1 due to changes in temperature 


and supply voltage, which affects frequency instability, will be discussed 
below. 


2 : é 

2 ) (9) — is the resolution factor for a cosinusoidal pulse with reference 
to a fixed component with parabolic approximation of the characteristics, 
which is necessary in designing the oscillator base circuits. 
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DESIGN RELATIONSHIPS IN STATIC CONDITIONS 


The design relationships for static conditions were obtained by the 
method of symbolic reduction of equations, developed by S.I. Yevtyanov [10]. 
In order that the effect of base current of a 
suitable oscillator (Fig. 1) is properly con- 
sidered in these reduced equations, the anal- 
ogous electron tube oscillator used for cal- 
culating grid current in reference [9] should 
be replaced by a generalized equivalent cir- 
cuit (Figure 2). 

The processes occurring in the circuit 
of Fig. 2, taking into account the accepted 
symbols, may be expressed by the following 
symbolic equations 


Figure 2 ; . 
U=Lan (:, re iy) e 


m 


(9) 


The impedances Ded and Z,, of the circuit in Fig. 2 are related to the 


impedances of Fig. 1 by the following expressions [11]. 


SE a TT Oe pea eA Ae) 
m . . ° ’ 0 5 = . 10) 
Z,+2Z2+Zqu 2,+ Z2+ Zqu ( 
Here 
2 = ig (latest) 
Tt 1+ ito(1 + by) (11) 


is the impedance of the crystal; 


the reactive part of the circuit (Figure 2). 
In Equation (11) and (12), the following symbols are used: 


Vi 20 5 Fe = 1] — is the generalized detuning of crystal relative to the 
Oy 
frequency of its series resonance a, [5], @, 4, and Q — are, respec- 


tively, the frequency of the generated oscillations, the series resonance 
frequency of the crystal, and the crystal Q - factor; r, — is the resistance 


of the dynamic part of the crystal [9]; To = ayry.Cg — is the crystal param- 


eter, taking into account the effect of its static capacitance Co. 
In expression (12) it was considered that the difference between the fre - 
quencies of self-oscillations «and of the quartz 4), may be neglected. 
Substituting (11) and (12) in (10), and then the latter result in (9), 
we obtain the following reduced equation showing the relationship be- 
tween the complex amplitudes of the first harmonics of the collector cur- 
rent Tg and base current Ip, with the amplitude of the base voltage 


Uo. 


ee | rg ee ey egg | Lee 
Ui Gna) OER fag aaetdealg So aati oe 


4+ 4+%K,)(1-+ i))}. (13) 


XX 
1 ; 
Here R = EE oS Se is the control resistance of the oscillator at a 
k 
frequency &}, without considering the capacitance Co; T Cte ryCy —is 
the parameter expressing the effect of the circuit-switching capacitance; 
12 : rong ed = 
Cr, = eee is the total capacitance of the circuit; Kg 
circuit capacitances. 

The expressions for calculating the amplitudes and frequencies of the 
oscillations under stationary conditions may be obtained from (13) if we 
adjust both parts of it separately into real and imaginary components. At 
the same time, it is necessary to take into account, by analogy with a tube 
oscillator [9], that 


— is the 


1 
Co 


a Uo (14) 


where S, and S,, — are the complex, quasi-linear transconductances of 
currents in the collector (6) and base (5). 

With the above considerations, the following expressions were ob- 
tained for calculating the frequency and amplitude of oscillations in a sta- 
tionary state: 


a Caw 9] B + {451 [¥ (1 — to») — to] + BY Sot (15) 
S [vy (1 — tov) — *0o] B — Fs i 
S 
S = = 16 
iRy(1 ls io ’ ( ) 
Bo 
where Soy, = Soil, ik, & — is the normalized base current curve; 
B=1+ at — is the coefficient which takes into account the base 
1 


current 


Ry=R{(1 — oe al 


(17) 
is the controlling resistance of the oscillator 
ire v(1 — 70 ¥) — to 
v2 4 (PK) Swi oe 
e hesakiene (18) 


— is the transmission factor, taking into account the effect of base current 
on the oscillation amplitude under stationary conditions. 

From (15) it may be seen that there are two values of generating fre- 
quency v for each value of circuit capacitance Cy, (or TL). The investiga- 
tion of stability in reference [12] shows that under "soft" operating conditions 
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(operation with self-bias), for each value of Cy, (T,) there will be created 
and sustained an oscillation at the lower of the Pee emencics v determined 
from Eq. (15). 

A detailed study of (15) and (16) also shows that in master oscillators 
using PPT (junction transistors), in which the collector current transcon- 
ductance characteristic S is the order of hundreds of ma/v, stable stationary 
conditions will be obtained with large circuit capacitances. In addition, the 
following condition is fulfilled: 


Sy Dy > yo (19) 


Besides, in the frequency band under examination, usually 


<1 (20) 


and, in the conditions recommended below, the following relationship is 
fulfilled 


i (21) 


Considering (19) -(21), we obtain from (15) -(17)the following approxi- 
mate expressions for calculating the parameters of circuits under stationary 
conditions: 


(22) 
GC => nN ee £ — a 
: aie V Ryrx 6 Gs Ko’ (23) 


where X — is the wavelength of the self-oscillations, G = a = ils 
the regenerative factor. Yq (2) 
With known parameters for the triode and crystal, the following order of 
' calculations of the oscillator is recommended: select a regenerative factor 
G > 1 and determine from (22) the required value R_. Next, a value should 
be assigned to K, and the capacitances C, and Cg should be calculated by 
formula (23). Then it is necessary to select the initial data for determining 
the operating conditions — collector voltage E), and the static component of 
the collector current I,,g. This method of calculation is similar to the cal- 
culations for the:tube oscillator in reference [9] and described in [6]. 
Information on a rational selection of G, Kg, E;, and I. are given below 
on the basis of investigations of frequency instability. 


METHODS FOR INVESTIGATING FREQUENCY INSTABILITY 


Frequency instability due to temperature variations is caused by 
changes in the parameters of the transistor, circuit, and crystal, 
whereas frequency instability due to variations in line voltage is caused 
by changes in the parameters of the transistor and by the power dis- 


sipated in the crystal. 
' From (15), the following expressions were obtained for calculating 


the components of frequency instability caused by changes in: 
a) output capacitance of the transistor 
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SAT ota fo ae) RL SOY V5 
t ( xe) 2 QxryC? a wel (24) 


b) input capacitance and conductance (base current transconductances) 


af es | ee i ene 
= (Coe) = — ————_ AC , = (S.)= Oly 
ip (Coe) 2 QxreC3 Me aS 2 QuC2 (25) 
c) phases of collector current transconductance 
- eke 
(P51) = cama ra (26) 

d) capacitances C, and Cy of the circuit 

BIN CY cee CE en ee ee 

a. pemm lars on ge 2QurxCo%e C2” (27) 


e) static capacitance and resistance of the dynamic branch of the crystal 


ef. = — eke se Seas 7 — Meh eCo_ Are 
f sis 2 QyreC2 wy Co’ ()= SER Is (28) 


The expression for calculating the instability components caused by 
changes in power dissipated by the crystal A Py y Was obtained by the method 
described in [5] and has the form 


Af (pe EK Faqu fo (81 OF 3 to 
a Kqu KoSi iP (29) 
where TKF ay, Ky yu — are the frequency temperature coefficient and the 
heat emission doetticient of anh crystal; according to data in [13] for non- 
vacuum type crystals ao 0-2 wt/deg; g1(9) — is the form factor for a 


cosinusoidal pulse with line- wees approximation of the characteristic. 
An analysis of the resulting expressions shows that the instability com- 


ponent ot (2g) does not depend on the capacitances Cy and Coy, and the 


components 


of of of oF of 
f (Cee) f (Coe), f (So), f (C,), f (C3) 


decrease as Cy and C2 increase. 

In agreement with (23), with other conditions equal, the capacitances 
Cj, and C, increase with decrease in R,,, that is, with a decrease in the 
regeneration factor G (22). In addition? when G decreases, the instability 


f 
factor ais (Pau also decreases (29) since the angle of collector current 
cutoff increases and, correspondingly, g1(8) drops and $1 builds up [9], 


; ny Ne 
while the ratio ae , as investigations show, hardly changes during 
ie 
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variations of G within wide limits. Therefore, during the calculations, it 

is desirable to select small values for G. However, in view of the dis- 

persion in transistor parameters and also possible fluctuations in line 

ae it is recommended that the selection of G be made in the range 
= 2-4, 

Recommendations on the selection of Kg, Ex, Iq may be worked out 
after a preliminary clarification of their effect on the frequency instability 
arising from changes in temperature and line voltage. 

The investigations of instability was carried out by calculation and ex- 
periment on an oscillator using a P-403 transistor (fy 100 Mc) at fre- 
quencies of f= 0.75, 2 and 5 Mc with operation in an undervoltage condition. 
Preliminary measurements showed that for an overvoltage condition the 
frequency instability rises sharply. 

The experimental oscillator was designed in accordance with the cir- 
cuit in Fig. 1 Provision was made for two methods of supplying power to 
the base: from a separate source (1) and from the collector voltage supply 


(2). 

The calculations of the instability components with changes in tempera- 
ture and supply voltage were carried out by the formulas obtained above. 
The resultant instability was calculated as the sum of the separate com- 
ponents and verified experimentally. 

The current Ij,g varies with changes in temperature and supply voltage. 
The increments in current AlIj,9 caused by changes in temperature or supply 
voltage, which are necessary in calculating instability, may be measured 
by the instrument I,o in the circuit of Fig. 1. 

In calculating the instability components under stationary operating 
conditions, it is necessary to determine ACje1, ACoe1, ASp,, AMs1. 
which depend on changes in temperature or supply voltage and the increments 
in current Al;,g created by them, with high amplitudes of voltage Ug and 
cutoff of collector current (changes in the quasi-linear parameters of the 
transistor). The approximate calculation method in [8] does not permit a 
determination of changes in the quasi-linear parameters for low-level os- 
cillations Ex, I, and t®°. Therefore, the changes in transistor parameters, 
which are essential in calculating instability, were determined from ex- 

- perimental relationships of Cye, Coe, So, gs (Ex, Ix, t°), obtained for 
Ug ~0 with the assumptions taken in reference [7]. These assumptions 
lead to the hypothesis that the relationships of Cre, Coe. Sg, Os (Ek, Ik, 

t°), obtained for Ug> 0 also hold true for high amplitudes if the current at 
the operating point with Uo > 0 is replaced by the fixed component of the 
collector current under stationary conditions. As in [7], experimental 
verification (Fig. 4, 5) showed that under these assumptions the calculated 
and experimental results conformed to a sufficient degree. 

Fig. 3 shows the experimental relations between Cye (Ek, Ik, t°), and 

“Coe; So, Pg Ug, t’) for the P-403 triode at f = 2 Mc. The dependence of 
Coe, So, Pg on E, was expressed somewhat weakly and, therefore, is not 
shown in Fig. 3. The measurements of the transistor parameters at fre- 
quencies of f = 0.75 Mc and f = 5 Mc show that the nature of the indicated 
relationships are maintained, although the absolute values of the parameters 
vary. 

vat all frequencies, with known variations in temperature or supply 
voltages, known increments of current AI,g and experimental relationships 
between Ce (Ex, Ix, t°) and Coe, Sg, Os (ent), Are, ACog, ASo,. Ag 
were determined and both the component as well as resultant instabilities 


were calculated. . bi 
During the experimental verifications the resultant frequency instability 
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was measured directly. 
In reference [7] it was shown that with no change in Ej, and Io and an 
increase in the self-bias resistor Re? in the emitter circuit(Figure 1), the 


Al;.g decrease due to a change in the temperature or supply voltage. Thus, 
with an increase in R,, the frequency instability due to supply voltage 
variations will decrease. The temperature instability of frequency in an 
oscillator with a drift transistor is determined basically by the changes in 
phase of the transconductance from temperature and the resulting incre- 
ments in current Al;,g [7]. To lessen this instability it is not expedient to 
select a much larger R,. The difficulty is that with an increase in tempera- 
ture, I,g starts to rise. In drift transistors, an increase in temperature 
and current Ij,g causes a change in phase of the transconductance of opposite 
sign (Fig. 3). Therefore, by selecting Re ~ Re opt it is possible to attain 
such increments of current Alj,g which will assure a reciprocal compensa- 
tion in frequency drift re- 
sulting from temperature 


af Ket K-05 4 variations in the triode 
parameters. 
a X, Ohms) In connection with the 


above, the effects of Re, 
as well as of Ko, Ex and 
Ig, on frequency insta- 


Er-4vY~ y= 905-10° 


: bility were studied. For 
-4-10 o = 
eaene: Oy Be this purpose, at each of 
the frequencies (0.75, 
Figure 4 2and5 Mc), during changes 


in temperature or supply 


>The voltages on the transistor collector and base differ from those of the 
voltage sources Eg and E (Figure 1) because of the voltage drops in the 
resistor R, from the collector current Ig and base current Ip9. There- 
fore, with a change in Re, the current Ig and E; will vary. In order to 
hold Ig and Ej, constant with increases in Re, it is necessary to raise the 
source voltages Eg and E accordingly. 
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voltage, the calculated and experimental dependence on Re was obtained of 


Af 


7 caused by changes in the transis- 


the resultant frequency instability 
tor parameters and power dissipated by the crystal. At the same time, the 
regenerative factor G was held unchanged; Ko, E, and Ig were varied. In 
addition, in order to be able to evaluate the effect of changes in power 
Af 

(Pqu) and 


f 
Re were obtained under these same conditions. The selection of Kop with 


dissipated in the crystal, calculated relationships between 


G = const was accomplished by varying the capacitances C, and Cy. For 
each set of values of Ko; E, and I.o> with a change in R,, the voltage 


sources Eo and E (Fig. 1) were simultaneously regulated in order to ob- 


tain, with G, Kg = const, the relation at (R,) with unvarying E. and 


Ixo- 
While investigating temperature instability, besides the indicated re- 
lationships, additional measurements were obtained of frequency instability 


caused by changes in temperature of circuit parts and crystal. 


FREQUENCY INSTABILITY DUE TO TEMPERATURE CHANGES 


At frequencies of f = 0.75, 2 and 5 Mc the variations in as (Pgu) 
Af 

if 
peratures of +30° C (from +23° C to +53° C) G= 4 for various values of 


and as Re varies were obtained with a variation in transistor tem- 


Ko, Ej, and Iyg. At all frequencies the component fe (Pqu) proved to be 


less than the resultant instability ( = ) and it may be considered that the 


latter is determined mainly by changes in the transistor parameters. 


As an example, Fig. 4 shows the calculated and experimental rela- 


. tionships of = (Re) obtained with f = 2 Mc, E;, = 4v, I, gq = 4 ma and 


various values of Ko: 


Fig. 4 shows a sufficiently close adherence of the calculated and ex- 
perimental relationships. This confirms the correctness of carrying out 
the calculations with the indicated assumptions. 

It may be seen from Fig. 4 that for all values of Ko it is possible to 

select Re = Re opt which results in the frequency being independent of 
temperature, i.e., a reciprocal compensation takes place of frequency 
~deift due to changes in the transistor parameters [7]. The resistance Re opt 
may be selected experimentally. 
The value of R, opt is only slightly dependent on Ko. When R, > Re opt 
the temperature instability of frequency is not large. 
Tests show that the values of Re ppt remain practically unchanged when 
E;, varies from 3 - 5 v, but are Hepemient on the current I,9. In varying 
Ig from 2 - 6 ma, a separate Re = Re opt could be selected for each setting. 
The investigations at f = 0.75 Mc and f = 5 Mc also confirmed the above 
results. 
’ The calculations and measurements show that with fluctuations in 
temperature of +1° C, the change in triode parameters (for R, > Re opt) 
cause frequency instability of the order of one part in 10-8, change in 
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capacitance of the circuit of 1079, change in capacitance C, and resistance 
r, of 10719, and the TKF gy of 10-6, These results coincide with data from 
experimental investigations in reference [3] and show that temperature in- 
stability of frequency is determined basically by the TKFgqy (the temperature 
coefficient of frequency of the crystal). In order to decrease it, the first 

step is to arrange thermostatic control of the crystal. Taking into considera- 
tion the small dimensions of the transistor and circuit parts, and the low 
level of power dissipated in them, it may be recommended that the whole 
oscillator be thermostatically controlled. 


FREQUENCY INSTABILITY DUE TO CHANGES IN SUPPLY VOLTAGE 


We have already investigated frequency instability during voltage 
variations on the collector E;,, the base Eg and source E in networks with 
common supply to the collector and base (Fig. 1, position 2). 

Investigations at frequencies of f = 0.75, 2 and 5 Mc showed that with 
G = 4, various E,, Ikg, Re > Re opt and Ky = 0.5 - 1.5, the change of 
one collector voltage by +20% causes a frequency deviation due to the 


There is considerably stronger dependence of frequency on voltage of 
the base Ep or of the source E, with a common supply, in which case the 
frequency instability from Eg and E are of the same order. This is ex- 
plained by the fact that with a change in E, ina common supply circuit, 
the voltages on the collector E), and base of the triode also change. The 
frequency stability is determined mainly by changes in the latter since its 
effect on frequency is greater than that of E,. Therefore, the results of 
the investigation of instability due to Eg relate in the same degree to the 
common supply circuit. 

Figure 5 shows the calculated and experimental dependence of At (Re) 


triode parameters of , and due to the power dissipated in the 


and the calculated dependence of — (Pqu) on Rg obtained for various Ko 


af 
hi 150 Ry (ohms) 
eh 4 


with f = 2 Mc, I,g =4 ma, Ex = -4v and a variation in Ep, (or E) of +20%. 
The figure shows a sufficiently close concurrence of the calculated and ex- 
perimental relationships and also, as in Fig. 4, confirms the correctness 
of the instability calculations on the basis of indicated assumptions. 

Figure 5 shows that frequency instability is determined basically by 
changes in the transistor parameters. Variations in power dissipated in 
the crystal are seen to be considerably less important. As Re increases, 


Af Af 
r and —F7 Pau) decrease. However, for large values of Reg 


(Re > 100 ohms), a further increase does not result in an appreciable 
gain in stability. Since the values of Re = Re opt are obtained on the basis 
of fairly broad conditions of temperature compensation (Fig. 4), then it 
is quite permissible to select Re > Re opt in order to increase stability 
during changes in Eg (or E), 


Figure 5 shows that instability of frequency due to Eg decreases with a 
rise in Ko: The reason for this is that, with other conditions equal, with 


; : F et msi Af 
an increase in Kg the instability components are (Cra) a ee (Poy) fall, 


and = (Coe) and (So) increase. Only the component —— (.) is 


f 
independentof Ko: Calculations show that with Kp > 1 the frequency in- 


stability is determined mainly by the components < (sg), ts (Coe); 


f 
AL (So) (At (Cie) and 2 


independent of Ko on the resultant frequency instability is determined by 


(Pau) are small). Since the first one is 


the advantages in the relationships between the components Af Oe) and 


At (So). Taking into consideration the signs in the expressions for cal- 


culating these components in (25) and (26), and also the fact that with an 


increase in Eg (or E) the current I,,9 rises, and in accordance with Fig. 3 


we have AC, > 0, 4S, > 0, we may conclude that ac ) and 
0e 0 if e 

Af 

f 

must exist a Ky = Ko opt at which the frequency instability will be a 


(So) cause frequency drift of opposite sign. Consequently, there 


minimum. As an illustration of the above, Fig. 6 shows the calculated 

relationships of at (Ko) obtained atf=2 Mc, G=4, Tea 4 ma, EK, =-4v, 
= 
values of R,. The figure 


varied +20% and two 


shows that for both values 
of R, the frequency in- 
stability was at a minimum 
with Ko opt = 8-10. The 
selection of Ko opt is not 
critical. With Ky) = 1 - 1.5, 
the frequency instability is 


not more than 20-30% that 
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obtained with Kg = Kgopt. An increase in Kg above Ky = 1 - 1.5 does not 


produce an appreciable gain in frequency stability. It is brought about 


oes : coed 
through the selection of G and 1,9, that is with a given Ur(Uo= 0) [6] ) 
to a decrease in amplitude of the collector voltage U, Sale aed where 
{K] an Sera is the modulus of the feedback coefficient. The oscillator 


V +i 
power developed by the transistor is also decreased. Therefore, it is quite 
satisfactory to select K, =1-1.5. <a 

With Rg > Reopt and Kg=1- 1.5, the frequency instability during 
changes in Eg (or E) of +20% is of the order of one in 10% 

Investigations of instability for various values of E},, Ig with f=2 Mc 
and at other frequencies (0.75 and 5 Mc) also confirmed the above results 
and showed that at these frequencies, also, it was permissible to select 
Re > Re opt and Koa iienls 5. In order to increase frequency stability it 
is desirable to supply power to the base circuit from a highly stable source. 

The investigations also showed that frequency instability due to Eo 
(or E) is nearly independent of the selection of E, and Ig if, for each 
value of E,, and Ig, we select Re > Re opt and K, = 1 - 1.5. Therefore, 
the recommendations on selection of Ek and I, values concur with those 
made in [7]. This means that for the P-403 triode, the selected values 
should be Ig = 2-6 ma and E, = 3-5 v. 


CONCLUSION 


A method was developed for the engineering design of a transistorized 
master oscillator with a crystal connected between the collector and base, 
for frequencies of f< 0.05 fy. This method may be used for the design of 
master oscillators with any type transistor. 

Certain design recommendations may be considered general in nature 
for all types of transistors. In order to attain high frequency stability in 
master oscillators of any type transistor, it is advisable to select under- 
voltage operating conditions, to have a regeneration factor G = 2-4, to 
arrange the self-bias from the emitter current (or combined) and to provide 
a highly stable power supply for the base circuit. Depending on the type of 
ae the recommendations may differ as to the selection of Ex, Ig, Ko 
and Re. 

For master oscillators using the P-403 transistor, the following values 
are considered suitable: Ex PSO V, lege coum: Ko = liam oland 
Re > Re opt. These same recommendations apply to the P-402 triode. 

It should be noted that for master oscillators using the diffusion-type 
triodes P-13 - P-15, in order to achieve high frequency stability it is de- 
sirable to select fn ho LOvy, I «9 = 5-8 ma and as large a value as possible 
for Rg. The latter requirement is based on the fact that it is impossible in 
these triodes to compensate for frequency drift when the triode parameters 
change with temperature [7], and in order to raise the temperature sta- 
bility it is desirable to decrease the increment in current Ixo Which is pro- 
duced by the change in temperature. It is also recommended orientationally 
that for these transistors K, ~ 0.3 - 0.5. For a more exact solution of 
the problem of selecting values of Kp, it is necessary to calculate the de- 
pendence of frequency instability on es with given values of E;, Ig and 
R,, and then to select a value of Kp which results in low instability. 

In conclusion, I wish to express my thanks to N.K. Bekasov for his 
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participation in the experiments. 


12. 


te 


REFERENCES 


Strutt Schaffenhauser. Experimentell und theoretische Untersuchung 
eines Frequenzstabilisirten Transistor Oszillator fur 8 Mggz. "Archiv 
der Elektrotechnischen Ubertragung, MALS IN@e. tithe 

Ludloff Awender. Uber Transistorszillatoren mit Quarzsteuerung. 
"Elektronische Rundschau,' 1958, No. 3. 

B.A. Timofeyev. Transistorized Oscillator with Increased Frequency 
Stability. Transistor Electronics in Instrument Production. Oborongiz 
1959. 

V.P. Zhukhovitskaya. On the Design and Comparison of Crystal- 
Controlled Master Oscillators, NDBS. 'Radiotekhnika i Elektronika,"' 
IN@o ils eGR» 

V.P. Zhukhovitskaya. Design and Comparison of Crystal-Controlled 
Master Oscillators. Candidates Dissertation, MEI, 1959. 

S.L. Kunina. Design of Low-Frequency Master Oscillators Using 
Junction Transistors. 'IVUZ — Radiotekhnika,'' No. 3, 1958. 

S.L. Kunina. Frequency Instability of Master Oscillators Using a 
P-403 Transistor. "Elektrosvaz','’ No. 7, 1960. 

B. Ye. Petrov. Equivalent Circuit of a Junction Transistor for High 
Sinusoidal Voltages. Semiconductor Instruments and Their Applica- 
tions (No. 4). "Sovetskoe Radio, '' 1960. 

S.I. Yevtyanov. Radio Transmitting Apparatus. Svyaz'izdat, 1950. 
S.I. Yevtyanov. The Relationship Between Symbolic and Reduced 
Equations. 'Radiotekhnika,'' No. 1, 1946. 

S.L. Kunina. Investigation of Crystal Oscillators Taking Grid Current 
into Account. ''Trudy MEI,' No. XXI, 'Radiotekhnika, '' 1956. 

S.I. Yevtyanov. Theory of Crystal-Controlled Master Oscillators. 
"Radiotekhnika,' No. 1, 5, 1949. 

M. Ye. Zhabotinskiy, P. Ye. Zil'berman. Relationship Between 
Frequency of Crystal Oscillators and Power Dissipated in the Crystal. 
"Radiotekhnika i Elektronika,'' No. 2, 1958. 


> 


Received May 10, 1960 


69 


ON THE FORMATION 
OF TELEPHONE DISTRICTS 


UNDER CONDITIONS OF UNEQUAL DENSITY 


I.M. Zhdanov 


In working out prospective plans for the development of city telephone 
networks, the layout of the telephone districts is carried out on the basis of 
the principle of equal-capacity districts. A brief explanation of this prob- 
lem is given and a new principle is proposed for the layout of telephone dis- 
trict s — the principle of equality of the total length of the subscribers' 
lines. By means of an example, the possibility of decreasing the cost of 
line construction is illustrated. Along with this, attention is also given to 
the necessity for exposing and calculating inequalities in district telephone 
densities, and through an example, the effect of inequalities on the length 
of the subscribers' lines, and on the degree of accuracy achieved in de- 
termining this length, is shown. 


INTRODUCTION 


In the prospective planning of district telephone networks many difficult 
tasks arise, the correct solutions of which may lead to a considerable saving 
in costs and materials, since the volume of work on district telephone net- 
works amounts to many millions of rubles. 

Among these tasks, by no means the least, concerns the determination 
of the correct principle for the layout of the telephone districts, and solving 
the problem of the effect of, and calculations for, the inequalities in dis- 
tribution densities of subscribers in the districts. The correct solution of 
these problems can bring about essential economies in resources. 

In recent times, during the preparation of the general layout for the 
telephone systems of Moscow (1953) and Leningrad (1960), and also in 
publications [1], the principle of equal-capacity districts (equal capacity 
for district ATS [automatic or dial telephone stations or exchanges]), as 
the basis for laying out telephone districts, has had considerable applica- 
tion. This principle, without having adequate basis, is put forth as the 
basic principle both for the equal distribution of subscribers (where its 
correctness is not disputed), as well as for unequal district densities, 
which can lead to an increase in material outlays, particularly where the 
variation in density is large. 
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As a rule, the subscriber densities in districts in the central part of 
a city's territory are higher than in the peripheral areas and reaches a 
ratio of 30:1, as can be seen from the general layout of the Leningrad 
telephone system. 

With such significant inequalities, the use of the principle of equal- 
capacity districts leads to a considerable rise in construction costs. 

The present article intends to throw light upon another principle for 
laying out telephone districts, with an without taking into account the in- 
equalities in density, and to give a comparative evaluation. Along with 
this, some results will be presented of preliminary investigations on the 
problem of the effect and calculation of inequalities. 

Everything, together with the preceding, will be directed toward maxi- 
mum clarification of the problem of district layout and to anticipate all 
possible essential deviations which may be permitted in laying out the 
telephone districts. 


PRINCIPLE FOR LAYING OUT TELEPHONE DISTRICTS 


As is known, the basic quantity in determining the number of districts, 
and for laying them out within the territory of a city of area S, is the optimal 
capacity of the district ATS — Cg. 

The number of stations or districts — n for an equal subscriber density 
within a district o (Figure 1) is determined from 


Co’ (1) 


where o — is the territorial density in T/ha [telephones per hectare]. 
In the particular case when 
S = 7200 ha, o = 25 T/ha and 
Co = 10, 000 subscribers, we have 


_ 1200-25. _ yg 


& 10, 000 
Yj The capacity of the total net- 
| work will be 
N =So = 7200-25 = 180, 000 
subscribers. 
In order to have the shortest 
total length L, of subscriber lines 


for all districts, the 18 district 
ATS should be distributed in accord- 
ance with Figure 2, which shows 
districts of equal area and capacity, 
with equal and minimal average 


values for the lengths of the sub- 
scriber lines — l;, that is, the 
following equalities are preserved 
2]: 

Vs 


oem) 


S,=S,=S,=...=S,, (2) 
CC 0 oe ee (3) 
Figure il Cy, = Cl, =C,/3;= oo (Ciel be (4) 


Cs 


where S;, Sg, Ss, etc. — are the areas of the territorial districts; Cy, Co, 
C., ete. — are the capacities of the districts; and /1, /5, lg, ete. — are 
the average lengths of the subscriber lines in the districts. 

Consequently, in laying out the telephone regions, we may be guided by 
the following principles corresponding to expressions (2), (3), and (4): 

a) equality of district territories, or equal-sized districts (2), 

b) equality of district capacities, or equal-capacity districts (3), 

c) equality of total length of subscriber lines in the districts (4). 

In the case of equal density, all three principles are equally correct. 

But in practice, we already have cases of unequal distribution of 
density. 

Figure 3 shows one such density distribution where the average density 
04, corresponds to the central area S, and the average value 09, to the 


Figure 2 Figure 3 


peripheral part with a total area So, in which case 0, > 09(Figure 3 shows 
these areas in different crosshatching). 


In reference [2] by this author, it was shown that in this case only one 
third of the principles is applicable — the equality of total length of sub- 


scriber lines. The remaining two principles will lead to a longer total length 
of subscriber lines L.. 


In preserving the third principle under the new conditions. E 4 
be transformed into the form ec Nagi 


SF, VS, =Sy20F2V_ Sq == Sy0,F V S3=Syo.F,V Sy, 6) 


where, in place of Cy; Co, Cy, etc. we have substituted S, o,, 
S992, SgGs, etc, and in place of Las lo, l., etc. we have taken 
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Pilesi. Ff; VS. F, V Sz ete. 


Herern Op wa Oi 0, — are the subscriber densities in the districts; 
S;, 5,, 5, — are the areas of the district territories, and F a Lionel ar 
etc. — are the coefficients which characterize the subscriber distribution 
in the districts, the location of the ATS in the districts, the shape of the 
territory and the direction of the main cabie layouts. The value of the co- 
efficient F for a uniform distribution of subscribers in the districts, with 
the ATS located in the center of a square district territory, and with the 
lines laid along two mutually-perpendicular directions, will be, as is 
known, 0.5 — and accordingly, 1 =0.5 Vs or [ = 0.5 a where a is one 


side of the square. The value of F for other rectangular territories, or 
approximations of them, may be found in the literature [3]. 

As can be seen, Equation (5) contains the value 0j which characterizes 
the densities by districts. 

Consequently, in using the third principle and its corresponding equa- 
tion (5), inequalities in district densities are taken into account, as a re- 
sult of which this principle will give a shorter total length of subscriber 
lines Lg. 

This is confirmed by the example below. 

Let us assume that 18 districts (n = 18) are laid out in a fixed territory 
S with uniform distribution of densities in accordance with Fig. 3 where 
the central area S;, has a density 0, which is greater than the density 05. 

Laying out the telephone districts on the basis of the first principle — 
equality of district territories (2), we will have a distribution in accordance 
with Figure 2. The equality of district capacity (3) gives a distribution as in 
Figure 4 and for the third principle — equality in the total lengths of sub- 
scriber lines — we have the distribution of Figure 5. 


n=10 ,=6 
n,=13 ny=5 Pees: 


Figure 4 
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Table 1 


Increase in length of 
subscribers' lines in 
comparison with the 
third principle 


Principle for laying out 
the telephone districts 


Principle of equal-sized district 
territories (Figure 2). 

Principle of equal-capacity dis- 
tricts (Figure 4). 

Principle of equality of total length 

of subscribers' lines (Figure 5) 


175, 440 


148,920 


138, 280 


1 5 er) : 
Data calculated with S = 12,000; 5S, = TED S; So = 3 83 Pe=0253 


Oras 48 t/ha; a Mi 4 t/ha. 


In Table 1, in all three cases the values for the total length of the sub- 
scribers lines L, were calculated by the formulas: 


a) for the 1st principle 


Lu = (Siar + Sie) F VS, 
: (6) 


b) for the 2nd and 3rd principles 


Fes ise 
1D ay ata V 5, + S,39F = R e 


The footnote indicates the initial data taken for calculations by the 
above formulas. 


The results of the calculations completely affirm the correctness of 
the expressed statement on the advantages of the 3rd principle of laying 
out telephone districts as the one which gives the shortest total length of 
subscribers' lines L.. 

The divergence in the values of L_ was essentially 27% for the 1st 
principle, 8% for the 2nd principle, which must be taken into account, 
particularly when it is considered that each percent of saving in cable ex- 
penditures for subscribers' lines may add up to a considerable sum. 

This divergence may be increased still further if the difference in the 
values for densities of the districts is increased. In our example, the 

ie pire Sieg eet ; ; 
a ; eer Ee 12, while in practice it may 
be higher than this. In regard to expenditures for station equipment, since 
the number of stations (districts) remained unchanged and the capacity of 
the entire network also remained constant, there is no basis for assuming 
that there will be any essential change in this factor for the three principles 
examined. 

From all of the above, we may draw the conclusion that for conditions 
of ‘unequal distribution of subscribers' densities by districts, when laying 
out the latter it is necessary to be guided by the third principle — equality 


density ratio was taken as 
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of total length of subscribers' lines. 

The principle of equal-capacity districts which was used in the general 
layouts for Moscow and Leningrad, as we can see from the above, does not 
give a correct solution of the problem. It inevitably leads to an increase in 
line expenditures, the amount of which will become larger as the divergence 
in the subscribers' densities in the districts increases. 

As regards the principle of equal-sized districts, being suitable only 
for conditions of equal distribution of subscribers with unequal densities in 
the districts, it may result in an even more substantial increase in cable 
expenditures for subscribers' lines than the 2nd principle. 


ON THE EFFECT AND CALCULATION OF INEQUALITIES 
IN DISTRIBUTION OF SUBSCRIBERS IN CITY TERRITORIES 


The inequalities in subscriber distribution in a city territory are 
two-fold: on the one hand, there is inequality of density (for average values) 
in the districts, which was examined above, and on the other, there is an 
unequal spacing of subscribers within the district. 

The first type of inequality, when using the 3rd principle for the Layout 
of telephone districts, is taken into account in (5) by the values of 01, 0.,, 
O., while the second type of inequality is also taken into acc aint in (5) 
through the values of F,, Boyt at rs 

For example, for a concentrated distribution of subscribers in a square 


territory (Fig. 7) it is necessary to take a lower value for a uniform dis- 
tribution (Fig. 6). 
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Figure 6 Figure 7 


The investigations carried out by P.P. Yushkov [4], N.Ye. [sic] Buk- 
hman [5] and also by the author of the present article [3, 6 ] (sic - 2, 3) shows 
that the change in the value of F arising from different distributions of sub 
scribers inside a district will be relatively small, for example, upto 50% in com- 
parison with F, either as an increase or decrease, will not result in any 
essential changes in total length of subscribers’ lines L,, inequality in sub- 
scriber distribution inside a district may be disregarded, i.e. a uniform 
distribution may be assumed. + 

Thus, it is necessary to take into account the inequality of densities 
of the districts which, as was shown in the previous example, may cause 
substantial changes in the value of Lg. ; 

In order to judge the relative importance of this inequality, another 
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example will be examined. 

Let us assume that a city territory with an area S$ is square-shaped 
(Fig. 8) and has N subscribers. Not knowing the actual distribution of sub- 
scribers within the territories of city S, we shall, for the sake of a simple 
and approximate solution of the problem, assume a uniform distribution with 
an average density value of Og, = 60 T/ha = 6000 T/sq km and have laid 
out 24 telephone districts (n = 24) in accordance with the pattern in Fig. Se 


Cav = 60 t/ha 
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Figure 8 Figure 9 
In such a distribution we have, calculating Ly by formula (6) 


La, = 6122S VY Skm. 


Further, for a more precise definition of subscriber distribution, 
let us consider that their densities in sections $1 and Sy are not the same 
(Fig. 10). Let us take, for example, a density of 0] = 10 t/ha for Sj and 
Og = 110 t/ha for So. 

Under such conditions the layout of the same 24 districts resulting 
from (5) should be carried out in accordance with Figure 11. 

Calculating the total length of the subscribers! lines by formula (7) 
gives 


Lyp=5258SV S. 


Continuing further a more precise and clearer definition of the in- 
equality, we shall assume, for example, that in section So (Fig. 10) there 
is also inequality in subscriber distribition in accordance with Fig. 12 


where together with the average value 05, we also have the different values 
Os, "ANGOigs 
34 


The calculation of L, in this case by formula 


acs . wd = 
Pes ME yj = + Sio, F y= Lisain Tse 
ny ng ng 


gives us a value of 


(8) 


L43= 4809S) S. 
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Further definition of this inequality (Figure 14) leads to a new distri- 


bution of these same 24 districts (Figure 15). Calculating La4 according 
to the formula 


stiaiity } pik acteurs Pars’ = 
Lay=S\9,F ye ote S} 0, F yt ale OF ok = +So,F y= 
Ty Ng Nz m 


gives, in this case, an even smaller value La4 = 477.28 Vs F 


(9) 


0,=10 t/ha 


S, ~__)=110 t/ha 


The above example shows graphically how the value of L, depends on 
the degree of definition of the inequality in subscriber density by districts. 

The more precisely the inequality is defined, the higher is the degree 
of accuracy in the calculations, and also the more correct can be the se- 
lection of the best possible variation for the layout and distribution of the 
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districts. Inourcase, the most suitable variant is the one shown in Figure 15. 
However, to simplify the problem and speed up its solution, the degree of 


0;=10t/ha 0,=78t/ha 
0} =66t/ha 5, ~_04=190t/ha 


ae 


Figure 14 


Figure 15 


precise definition of the inequality may be limited if further definition does 
not result in a substantial change in La. 

In our example it is sufficient to take three gradations of density 
(71, To, ), without resorting to further definition of subscriber distri- 
bution. Suck a limitation in accuracy of calculations will give an error for 
L, not greater than 10%. 

A braoder and deeper investigation of this problem by the author showed 
that the same limitation of three gradations may be applied to many other 
cases. When this solution must be carried out with greater accuracy, it is 
necessary to increase the number of density gradations K to the number of 
districts n, and only with k = n can we be assured of a minimum L.. 

The above examples enable us to set forth the following. Resorting to 
the most rational principle for laying out the telephone districts — the 
principle of equality of the total lengths of the subscribers' lines when 
laying out and organizing the telephone districts for a specific network ca- 
pacity and a specific number of district ATS — we shall obtain a total 
length of subscribers' lines L, which will be shorter for unequal densities 
in the districts than in the case of a uniform distribution. This can be seen 
graphically in the distributions of Figures 9 and 15. 
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This difference will be quite substantial when the subscribers are con- 
centrated in the central part of the city territory, in which case, aside from 
a decrease in total length of subscribers' lines, there will also be a decrease 
in connecting lines due to an increase in density of the ATS spa cing in the 
central part. 

Accordingly, inequality in density distribution in the districts creates 
more suitable conditions for a district layout of the telephone network than 
does a uniform distribution or equal values for the densities. However, this 
advantage may be attained only by observing the following conditions: 


a) equality in total lengths of the subscribers' lines for the laid-out 
telephone districts, 

b) sufficiently-detailed definition of the inequality in densities by dis- 
tricts, or by larger areas containing several districts. 


In practice, the calculation of the inequality and laying out the tele- 
phone districts in accordance with established norms of line construction 
costs, for subscribers' lines, may be carried out in the same manner as is 
done in the general layouts ''Genskhema" except that the principle of equal- 
ity of total lengths of subscribers' lines is used instead of the principle of 
equal-capacity districts. 


CONC LUSION 


The limitations of this article do not permit a more detailed examina- 
tion of this problem and thus, to deal more thoroughly with the theoretical 
aspects. The latter was covered to some extent by the author in reference 
[2], to which particular attention is called. 

But even in the present article it is obvious how important it is to 
have the right approach to the solution of the problem of laying out telephone 
districts, taking into account the inequality of subscriber distribution. 

On the basis of the examples given and conclusions reached for each 
particular problem, we come to the following general conclusion. 

The principle of equality of total lengths of subscriber lines is the best 
one, and with fixed costs for station equipment, will result in minimum ex- 
- penditures for subscriber-line cables. 

The opposing principle used in practice, that of equal-capacity districts, 
can compete with the proposed principle only in the case of equal densities 
in the districts, and even then is permissible only for very insignificant de- 
viations in densities. 

With density ratios of 1:10 and higher, the presently-accepted principle, 
as was shown in the example, will lead to significant increases om cable 
_ expenditures of the order of 8% and higher, depending on the degree of de- 
viation in district densities. 

Nonuniform distribution of densities, particularly with concentration 
of density in the central part of a city territory, creates the most favorable 
conditions for laying out of a telephone system, more so than a uniform 
distribution. 

This advantage may be realized with sufficient nonuniformity in the 
district densities. The above situation should prevail in each district sepa- 
rately. In some cases, when a high degree of accuracy in the determination 
of L,, is not required, an average value of density may be used for a group 
of districts, keeping in mind that the approximation of the most desirable 
solution depends on the density gradations, the number of which must not 
be less than three. (With three gradations, the error will not exceed ap- 


proximately 10%). 
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The basic quantity in laying out districts, as is known, is the optimum 
capacity, on the basis of which the number of districts is determined, and 
which also determines the layout for the telephone districts in accordance 
with the principle used for laying out the districts. The optimum value for 
capacity with different densities in the districts need not remain fixed, and 
will increase as the district densities increase, and vice versa. 

The principle of equal-capacity districts makes possible an equal- 
capacity arrangement of district ATS, which is occasionally pointed to as 
an advantage of this method. 

From our viewpoint, such an advantage (if it may be considered an ad- 
vantage at all does not sufficiently compensate for the increase in cable 
costs which of necessity results in comparison with the recommended 
principle of equality in total lengths of subscribers' lines. 
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TRANSMISSION PARAMETERS 
AND EFFECTS OF TYPE TO CABLES 


IN A FREQUENCY RANGE TO 600 KC 


A.D. Ruga and Ye. V. Kostyanets 


The transmission parameters and effects of using Type TG cables in 
a frequency spectrum to 600 ke are presented. Conclusions are reached 
regarding the possibility and advisability of hf (high frequency) multi- 
plexing of circuits in existing interoffice connecting cab es of GTS [City 
Telephone Networks]. 


INTRODUCTION 


The intense development of city telephone networks (GTS) is tied in 
with considerable capital investment and large outlays of such short-supply 
metals as lead and copper. One of the most effective methods for utilizing 
city cables is that of high-frequency multiplexing of the trunk-line circuits. 
; At the present time the GTS trunk lines are designed for MKS cables 

with multiplexing of their circuits by means of type KRR-30/60 equipment. 

A large, as yet unused, reserve is available through the multiplexing of the 
existing type TG cables of GTS. Investigations carried out in 1957-1958 by 
the TsNIIS [Central Scientific Research Institute of Communications ?] on type 
TG cable in the high-frequency spectrum, and also the investigations of 
NIITS [Scientific Research Institute of Telephone Networks ?] in recent years, 
show the feasibility of multiplexing in circuits using the existing type TG 
- cables, and data is also available abroad regarding the development of hf 
equipment for multiplexing GTS cables [1, 2, 3] which confirms the fact 
that sufficiently-low cost equipment may be produced for this purpose. 

The present article contains basic results achieved by NIITS of the 
Ministry of Communications in their investigations on the parameters of 
GTS cables of type TG in the frequency range to 600 ke. 


TRANSMISSION PARAMETERS AND EFFECTS OF TG CABLES 


In city trunk lines TG cables of various diameters and differently - 


designed conductor insulation are in use. 
Formerly, the industry produced, for the most part, cables with a 
paper-air insulation, type ''Trubochka"' [little tube], chiefly in diameters 
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0. 
0. 


but also 0.6, 0.7 and 0.8 mm. Cables with a conductor diameter of 


2 
ss . ' . 
4 mm were produced with insulation of type "'Treugol nik" [triangle ]. 


At the present time, cables of 0.56 mm conductor diameter are pro- 


duced with paper-composition insulation. 


x \d=0.5mm 
d=0.5mm paper composition 


Figure 1 


Figures 1, 2 and 3 show data on the‘transmission parameters of cir- 
cuits using these cables in factory lengths which represent the average 
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values of the index of 
characteristic impedance 
[Z,,] and its angle @,, and 
also the attenuation con- 
stant 8 and phase-shift 
constant @. 

An analysis of meas- 
ured data showed that the 
capacitance of circuits 
with paper-composition in- 
sulation exceeds that of 
circuits with type"Tru- 
bochka" insulation by 
11.2%. This explains their 
relatively high kilometer 
attenuation (on the average, 
0.1 - 0.2 nepers/km). 

Separate circuits using 
the same cable in different 
factory lengths vary some- 
what in attenuation, which 


Table 1 Tablee 


Frequency 500 ke 


Effective 


capacitance Type of 


insulation 


Minimum 
Average 
Maximum 


Paper 
composition 


Minimum 
Paper-air 


Average 
Beef ; (Trubochka) 


is explained basically by 
deviations in their effec- 
tive capacitances. The 
latter, for cables with a 
conductor diameter of 
0.5 mm, are illustrated 
in Table 1. 

The suitability of a 
cable circuit for multi- 
plexing is determined to 
a considerable degree by 
the extent of shielding at 
the far end (b,) and by the 
near-end crosstalk (Bj). 

Investigations of the 
dependence of b, and Bg 
in factory lengths, on 
frequency showed that the 
average value! of Bo in 
Figure 3 the frequency range from 

20-200 ke drops approxi- 
mately 2.5 nepers and with 
further increase in fre- 
- quency to 600 ke, remains nearly constant at a level of approximately 9 
~nepers (Fig. 4). The dependence of the average value b, on frequency is 
similar, decreasing by about 2-2.5 nepers in the frequency range to 300 ke, 
and then remains practically constant at about 8 nepers with further in- 
erease in frequency to 600 ke (Fig. 5). 
Figure 6 shows a cross-section of a TG-50x 2x 0.5 cable. As can 
be seen from Fig. 5 and 6, the circuit shielding increases with increasing 
separation of the pairs from each other. The GTS trunk lines, as a rule, 


1The values of By and bz (Figures 4 and5) were obtained as the result of averag- 
ing 33 identical pair combinations measured on several factory lengths. 
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Figure 4 


4.0 ee average value of shielding 
tO control wire-22 mF {3 shielding of separate pairs 


subjected to le---- {0 
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Figure 5 


are of the multiple type 
and contain conductors of 
different diameters, which 
make it expedient to char- 
acterize them by the oper- 
ating (effective) parameters 
by and Zin: 

The investigations were 
carried out on trunk lines 
extending from 1.7 to 10 
km. An analysis of the 
results obtained in meas- 
uring the effective attenu- 
ation shows that the devia- 
tion in attenuation, as with 
cables in factory lengths, 
is caused by deviation in 
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the effective capacitance, and does not 
exceed 0.25 - 0.3 nepers at 600 kc. The 
average value of effective attenuation is 
practically the same as the character- 
istic attenuation (Fig. 2). The dependence 
0 40 80 120 160 200 %0 260 320 of the scalar input impedance on fre- 
quency is of a pulsating nature, and at 
Figure 9 600 ke the deviation in the value of [Z;,] 
does not exceed + 15 ohms. 

An example of the dependence of 
far-end circuit shielding on frequency is shown in Fig. 7 where the average 
value is the arithmetical mean of all measured pair combinations in the line, 
and the minimum and maximum values are the arithmetical means of the 
measured combinations corresponding to the worst and best cases, respec- 
tively, in regard to mutual-shielding of the pairs. 

Figure 8 shows curves of average values of circuit shielding at different 
frequencies. 

Figure 9 shows the dependence of near-end crosstalk attentuation on fre- 
quency between two inter-exchange trunk-line cables 10 km in length. As can 
be seen in Fig. 9, the average value of crosstalk attenuation at 300 kc is 14 
nepers, which is ample proof of the practicability of using a two-cable system. 


DETERMINATION OF NUMBER OF PAIRS SUITABLE FOR 
MULTIPLEXING AND OPTIMUM MULTIPLEXING SPECTRUM 


, An analytical prediction of the number of pairs suitable for multi- 
plexing in low-frequency cables is usually carried out as follows [5]. 

With a total number of pairs in a cable of n, the number of two-pair 
combinations is equal to 


N = 1/2 n (n — 1) (1) 
If we assume that the number of pairs satisfying the given norms for 
shielding is equal to k, then, correspondingly, the number of possible two- 
pair combinations of multiplexed pairs analagous to (1) will be 


K = 1/72 k(k — 1) (2) 


The total number M of two-pair combinations in an n pair cable which 
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satisfy the shielding norms (blank squares in Fig. 10) will satisfy the re- 
lation ship 


Be. 
n 


M 


which is an indication of the low probability that the percent of suitable 

pairs will be equal to the percent of suitable combinations. These M - 
combinations may be arranged in the diagram (Fig. 10) by different methods, 
the number of which will be 


Seah CATT (4) 


where M may have any value from 0 toN. 
Assuming that the appearance of all arrangements of favorable combina- 
tions are of equal probability, the probability of the appearance of each such 


picture, which from now on we shall call constellations, may be determined 
from 


Joa MIN=M ye) 
Cc M ‘ (5) 
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Let C represent a single [individual] constellation, ay, — the number 
of combinations of k pairs in the given constellation, A, — average number 
of combinations of k multiplexed pairs, then 

PAE ee 
ry ( M ) ( 


M 


A= DY o.v(c)\=*e—” >) a. 


c=1 ic=1 


(6) 


The sum of the terms of the right hand member — which is the sum of 
all the suitable combinations of k pairs for all possible constellations — is 
calculated in the following manner: an individual k pair combination will be 
s uitable if all of its K combinations lie within the M category, which satisfy 
the shielding norms independently of the distribution of the remaining 
(M — K) combinations satisfying the norms. These remaining suitable 
(M — K) combinations may be arranged in various ways (constellations) on 
the diagram and their number will be 


Neh ee NEA 
M—K (M — K)! [(N— K)|—(M— K)\] © (7) 


The single k pair combination which we are considering will, therefore, 

: . aN —K 

be found in Cir igs 
ponents in the summation of Eq. (6). 

On the assumption of equal probability of the appearance of all con- 

stellations, all combinations of k pairs will also be of equal probability 


constellations and this number is the number of com- 


ee ley 
< kl(n—k)! * (8) 


Thus, we finally have 
Be (lil Hey hh): (9) 


If Mand N > K (the case of high frequencies), then 


Nak = Me K 
ye op rea bac eee 
In this case, expression (9) may be replaced by the approximation 
n\ ok 
aia (3)? (10) 


where p = oe 100 — is the percent of combinations satisfying the 


shielding norms. ; : 
On the basis of (10), a nomograph (Fig. 11) was obtained which permits 


the determination of the number of k pairs suitable for multiplexing as a 
function of the percent of suitable pair combinations p in the cable for 
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various values of n. 

It is not difficult to persuade oneself that the value k depends very 
little on the selection of the value B;,; therefore, the nomograph is calcu- 
lated for By = 1. The values of k obtained from the nomograph should be 
rounded out to the nearest, lowest whole number. 

Thus, it is obvious that in determining the number of pairs suitable 
for multiplexing, it is necessary to establish, experimentally, the value 
p, which will be different at different frequencies, types of cables and 
other factors. 

As was shown by comparing the experimental and calculated data, the 
error in an analytical determination of the number of pairs is small in 
practice, and may become significant (up to 50%) only in the case of cables 
with a small number of pairs (20-30). In addition, the analytical method 
does not exclude the necessity of measurements because if, for any par- 
ticular cable, p is known and k is calculated, then it is still necessary to 
carry out measurements for determining the actual pairs. 
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For the same cables in which the transmission parameters and effects 
were investigated, the number of pairs suitable for multiplexing was also 
determined. 

This determination was made on cables of different capacitance and 
length, and in each case a 50 pair group was selected arbitrarily for the 
investigations. The mutual shielding was determined for each pair with 
respect to each of the others. 

Measurements were carried out at frequencie 
dogs bones Sane quencies of 60, 100, 200, 300, 

As a result, curves were obtained of the dependence of the average- 
expected percent of pairs suitable for multiplexing on frequency and shielding 
norms for a repeater section of the line (/ = 3.5 —4km, TG-05 cable) 
The average-expected values of percent of pairs is shown in Fig. 12. , 

The determination of the number of the circuits to be subjected to 
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multiplexing is done by selecting them on the basis of the measurements 
of mutual shielding of each circuit with respect to all the others: 


N 100 
M FS 
% 
a) in the case when the average-expected number of multiplexed cir- 
cuits determined from Fig. 12 is substantially above requirements, it is 


necessary to measure the mutual shielding of M circuits, where N is the 
number of physical circuits required for multiplexing and the % refers to 
the average-expected percent of pairs (Fig. 12). These M circuits should 
be selected with maximum separation between each other within the cross- 
“section of the cable; 

b) in the case when the required number of physical circuits ap- 
proaches the average-expected value, it is necessary to measure mutual 
shielding of all circuits in the cable. 

In order to determine the effect 
of the capacitance of the measured 
group on the general conclusions 
regarding the number of pairs suitable 
for multiplexing, tests were conducted 
on these pairs for groups of different 
capacitances from 30 to 100 at 600 kc, 
which showed that the percent of pairs 
suitable for multiplexing does not de- 
pend on the capacitance of the group. 

Increasing the usable frequency 
range, while increasing the number 
of channels K (Fig. 13) obtained as 


0 00 20 300 400 500 Fc a result of the multiplexing of each 
circuit, causes, at the same time, a 
Figure 13 decrease in the number of circuits 
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suitable for multiplexing. In order to express the effect of these factors on 
the total number of channels obtained with the multiplexing of all suitable 
circuits in a 50 pair group, the curves in Fig. 13 were constructed as the 
result of measurements carried out on six 50 pair groups. These curves 
were constructed with the assumption that each channel occupies a band- 
width of 16.5 kc, which corresponds approximately to the frequency allo- 
cation in the KRR 30/60 equipment. 

Figure 13 gives an indication of the feasible multiplexing spectrum in 
relationship to the shielding norms, and shows the dependence of the num- 
ber of channels on frequency and length of line for a 50 pair group. 


CONCLUSIONS 


The results of the measurements carried out on factory lengths of 
type TG cables, and on interoffice trunk lines, showed that: 

1. There is considerable deviation in the values of the transmission 
parameters, which are caused mainly by variations in the values of the 
effective capacitances. 

2. The dependence of the transmission parameters and effects on 
frequency is of a non-monotonous nature, which can be explained by the 
presence of a large number of heterogeneities along the length of the line. 

3. The magnitude of near-end crosstalk between two cables at f = 300 ke 
is 13-14 nepers, on the average, which is evidence of the feasibility of 
using a two-cable high frequency system. 

4. Type TG cables contain a certain number of pairs suitable for 
multiplexing. Thus, for example, at 300 ke and with a shielding require- 
ment of b, = 7.5 nepers, it will average 27%, but at 600 ke it drops to 
7%, and in practice these values are independent of the capacitance of the 
group in which the suitable pairs are located. 

5. At existing prices of the KRR 30/60 equipment, its use in GTS 
[City Telephone Networks]? is quite suitable at distances of approximately 
12 km [4]. However, taking into consideration the substantial economy in 
non-ferrous metals resulting from multiplexing, the use of the KRR 30/60 
equipment will be justified in a number of cases at shorter distances. 

6. In cases where the KRR 30/60 equipment is used for high-frequency 
multiplexing, the approximate lengths of the repeater sections are shown 
in Table 2. 


Table 2 


Diameter of 
current-carrying |dq = 0.4mm 
conductors 


paper- | paper-air 
composition| insulation 
insulation Trubochka 


ae mm |dg = 0.5 mm 


*The use of KRR 30/60 equipment with type TG cables is possible only if 
the compensating network is changed accordingly. 
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